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Effects of turbulence on the received wave are studied using the 
Kolmogorov, Von Karman, and Gaussian models, and compared. The Gaussian 
model is simply related to the other two for purposes of comparison and 
is found to be a useful engineering approximation. An equivalent struc- 
ture constant is also defined for the latter. 

A generalized theory of the aperture antenna in turbulence is 
formulated. The uniform and tapered circular aperture antennas are 
specifically evaluated for any arbitrary angle of arrival of the wave 
in the main beam, and related to each other. 

The receiver, which has hitherto been neglected, is shown to play 
a fundamental role in the interpretation of signals which have passed 
through turbulence. Receivers are classed as synchronous and asynchron- 
ous. The signal measured by each class is described. Properties of the 
turbulence may be directly related to the output of the synchronous 
receiver defined in this study. 

Using experimental results, an equivalent atmosphere, homogeneously 
turbulent over the propagation path, with appropriate parameters, is 
shown to be suitable for prediction of long term average behavior. It is 
described by a height 6 km over an earth of 4/3 the physical radius, or 

8479 km. The long term average refractive index variance is 0.4x10“^^. 
The average correlation length is 46 m. The average equivalent struc- 
ture constant is 5.3x10"^'^. Extreme conditions are also estimated. The 
model predicted the experimental results over a range of frequencies of 
15 to 1. 

The variance and the mean signal level degradation as measured by 
the two receivers, the fluctuating power, antenna gain, and gain degrada- 
tion, are studied as functions of frequency from 3 to 90 GHz, aperture 
size, and elevation angle or equivalent path length. While the numerical 
results should be used with caution, the effects of turbulence are seen 
to be significant to the systems designer at the higher frequencies, 
lower elevation angles, or larger antenna apertures. 

Asynchronous receivers may be preferable for communications links. 
However, properties of turbulence are better measured by links employing 
synchronous receivers. 
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CHAPTER I 


INTRODUCTION 


1.1 Overview 


The allocation and use of ever higher frequencies in the microwave 
and millimeter wave regions of the electromagnetic spectrum has been 
spurred by the information explosion that followed the dawn of the space 
age. The attendant satellite communications technology, which had its 
modest beginnings in the ECHO satellite of 1960, generated its own 
demand. Taking over a large share of intercontinental telephone communi- 
cations, it also opened the way for wholly new services such as computer 
links, data nets, and intercontinental television links as well as 
direct television broadcasting by satellite. 

The resulting congestion has mandated the use of higher frequen- 
cies in hitherto unexplored regions of the electromagnetic spectrum. The 
new bands have been primarily determined by the characteristics of the 
terrestrial atmosphere. They are usually situated in 'electromagnetic 
windows', i.e., regions of relatively low attenuation away from resonance 
regions of its constituents, except in special cases such as secure 
short hop links. The fundamental constraint imposed by the windows on 
the choice of frequencies for general purpose usage precluded incremental 
changes in frequencies used and a gradual evolution of the technology. 

The new bands usually require quantum jumps in frequency along with a 
wholly new range of generating and receiving equipment. 
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The initially low powers and sensitivities, coupled with the usual- 
ly greater interaction with the atmosphere pose many problems to the 
systems designer. The effects of scintillation caused by atmospheric 
turbulence and bulk effects such as layers and fronts may all take their 
toll and need to be quantified. Attenuation statistics in clear air and 
rain are needed to ensure reliable links and to optimize parameters. 

One 'easy' solution to the problem of ensuring adequate system 
margins has been the use of large antennas, where permissible, to in- 
crease the directive gain. This increases the effective sensitivity of 
the receiver. However, it is now realized that perturbations of the 
received wave front across a large antenna by turbulence will, in effect, 
decrease its gain. This may also be thought of as an aperture to medium 
coupling loss. Indeed, at some point, an array of smaller antennas 
feeding coherently summing receivers begins to provide superior perform- 
ance at lower cost over a single receiver fed by a large antenna. 

The effects of the atmosphere on receiving systems had been 
encountered much earlier by optical astronomers. The familiar twinkling 
of the stars due to atmospheric turbulence, while celebrated in fiction, 
results in the degradation of the image at the focal plane of a tele- 
scope. The image spot may move randomly in the viewing field, causing 
a blurring of the image. In some cases, the image may even break up 
into several spots. These phenomena limit the resolution of the instru- 
ment. Further, the time available with good viewing conditions is also 
curtai led. 

Another related problem is the study of the surface tolerance of 
large antennas. Deviations from the optimum geometry of the surface due 
to tolerances in construction, thermal expansion, and the distortion 
caused by gravity result in a reduction in the gain of the antenna. The 
distortion of the incident wavefront by these deviations, when viewed by 
the feed system, is an exact analog of turbulence induced distortion and 
may indeed be indistinguishable from it. 
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Unfortunately, much of the study that has ensued is not in a form 
that is directly usable by the microwave engineer. The links between 
the parameters used by the optical astronomer, the microwave engineer, 
the antenna engineer, and the specialist in propagation have rarely been 
made explicit. Arduous measurements made by each group have usually 
remained within their own domains. The system designer has consequently 
relied on rules of thumb, empirical factors, overdesign, or has ignored 
the problem altogether. 


1.2 Outline of the Study 

References [11, 36, 44] have examined some aspects of the problems 
discussed above. The book by Shifrin [11] has been particularly useful. 
It provides a detailed analysis of linear and rectangular aperture an- 
tennas in the presence of turbulence. Unfortunately, there is very 
little material on the physical properties of the turbulence itself and 
the emphasis appears to be on the troposcatter problem. Theobold and 
Hodge [36] proposed a phenomenological model of the communications link. 
The wave incident on the antenna is decomposed into an amplitude varying 
component which is fixed in direction and a component which is fixed in 
amplitude, but varying in direction. Knepp [44] discussed a Gaussian 
tapered antenna. However, the results do not appear to be readily usable 
by the engineer. The last two works discuss on-axis incidence of the 
incoming wave only. None of the studies recognize the effect of the 
receiver on the measurement of the received signal statistics. 

Much of the literature referred to in this work are of use in a 
particular section only. Therefore, they are reviewed in the beginning 
of the relevant chapter and will not be presented at this point. 

It was felt, in light of the preceding discussion, that a coherent 
solution to the problem of a receiving system in the earth's turbulent 
atmosphere would be useful. Elements of propagation theory, of use to 
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the designer, have been selected from the wealth of information avail- 
able. A statistical analysis of an antenna in a turbulent medium is 
presented. The effects of the receiver on the signal gathered by the 
antenna, which appears to have been ignored hitherto, were analyzed 
with surprising and far reaching implications. The results are brought 
together to form a complete picture of the behavior of the communica- 
tions link in turbulence. 

This study attempts, first, to formulate the basics of the problem 
of the receiving aperture in a homogeneous turbulent medium. The medium 
is described using a simple statistical model. The concepts of element- 
ary antenna theory are used. The effects of the feed system are 
included as an integral part of the analysis. 

The propagation medium is then analyzed and several models of at- 
mospheric turbulence are compared. The statistics of the wave received 
at the antenna after propagating through the atmosphere are formulated. 

An analytically convenient model is then chosen and shown to be a reason- 
able engineering approximation to more commonly accepted models. 

The role of the receiver in the perceived signal statistics is 
studied next. The type of demodulator circuit used is shown to play a 
vital part in, measuring the fluctuations of the received signal. Thus, 
two general classes of receivers, called synchronous and asynchronous 
receivers, are defined. The statistics of the signal as measured by 
these receivers are formulated. 

The particular case of the circular aperture is then examined in 
greater detail. Uniform and Gaussian shaped feed illumination are con- 
sidered. The former also serves to link the results to the optical tele- 
scope in the particular case of focusing systems. The general behavior 
of the gain of the circular antenna in turbulence is shown. 

Next, experimental results are used to ascribe numerical values to 
the parameters of the propagation model. This describes the long term 
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average behavior of the atmospheric link. Estimates of the extreme 
conditions of turbulence are also made. 

Using these estimates of the average and extremes of atmospheric 
turbulence, the behavior of the variance and the average non-fluctuating 
power of the signal measured by asynchronous and synchronous receivers 
is studied. The degradation of the gain of the receiving antenna is 
also examined to obtain a deeper understanding of the effects of turbu- 
lence on microwave and millimeter wave satellite communications systems. 
Finally, some useful design curves based on these results are given, 
with extensions to some frequencies of future interest. 
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CHAPTER II 


THE ANTENNA 


This chapter examines the basic theory of the aperture antenna as 
a transmitter or receiver in the presence of statistically fluctuating 
media. The effect of the feed pattern is included. 

The far field pattern of an antenna acting as a transmitter is 
discussed briefly. The terminal voltage of a receiving antenna is then 
evaluated, first for plane wave incidence, to relate it to the case of 
the transmitter. The theory is then extended by considering the inci- 
dent wave to be perturbed arbitrarily in amplitude and phase in space. 

Finally, the perturbations are allowed to vary randomly in time 
and space to represent the effects of turbulence. The results are also 
related to the field in the focal plane of a focusing antenna. This 
links the antenna and the optical telescope. 

The aperture of the antenna is assumed to lie in the x-y plane and 
is denoted by S, as shown in Figure 2.1. P = (x,y) or P(IT), where R is 
a vector on the surface x-y, represents any point on the aperture, n is 
a unit vector normal to S. Q(r) is any point in the forward half-space 
at which the transmitted field E^ or the received field is measured. 
The time convention is used in describing phasors. 
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2.1 The Antenna as a Transmitter 

In this section, expressions are derived for the far field pattern 
of an aperture antenna acting as a transmitter for comparison later with 
the receiving antenna. 

Let the antenna produce a field distribution on the infinite x-y 
plane. Using the surface equivalence principle the aperture fields 
could be replaced by electric and magnetic equivalent current density 
distributions Jg(x,y) and fTg(x,y), where 
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Jg = n X IT 

Mg = E X n . 

f and TT are the electric and magnetic fields in the aperture plane. 
Alternatively, an equivalent electric current density 'J^{x,y) alone may 
be assumed such that 

= 2 n X IT . 

The above formulations are valid only in the half-space faced by the 
antenna. 

Further, let 

= Jq Pf (2.1.1) 

where Jq is a complex constant and f(x,y) is the dimensionless normal- 
ized aperture illumination function or feed function, such that 

and p^ is the unit polarization vector of the aperture illumination. 

Then, the radiated electric field pattern in the forward direction 
in the far (Fraunhofer) zone is given by [1], 


{QJ) 


“ j kr* 

f — jksine(xcos4) + ysincj)) 

2xr J e 


dxdy 


(2.1.3) 

where r,0,^ are the spherical coordinates of the observation point. 
is the characteristic impedance of free space. A is the wavelength, and 
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k, the propagation constant, is given in magnitude by 



The signed constant, k, is also the magnitude of T< where 


= k sine cos<t> 


ky = k sine sin<p 
k^ = k cose . 


Then the component of k lying in the aperture plane is 




(2.1.4) 


and 


Jno6 


-jkr 


E^(l<x,ky) - 2xr 


r T(x, 


y) e 


j(kj^x + k^y) 


dxdy. (2.1.5) 


Let 


g(l<x.V = f T(x,y) 


j(kxX + k^y) 


dxdy. 


(2.1.6) 


We recognize that g(k ,k^,) is the Fourier transform of the aperture 

^ y 

current distribution function. It is called the (unnormalized) electric 
field pattern function of the antenna in free space. Therefore, 


jn e 

EtC'x-V ' “if? '^0 9<''x'V • 

The magnitude of the field in any direction of polarization, p, in 
the far zone is given by 
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(2.1.8) 


JHoe 


•jkr 


~ 2Ar 


'’oP 


gO<X.I<y) 


Most of the transmitted energy is contained in the main beam of the 
antenna and the side lobe levels are assumed to be low. Consequently, 
the fields in the x-y plane and the equivalent surface density, 3^, may 
be considered to be small except over the physical aperture, S, itself. 
Therefore, the integration may be carried out over S alone with small 
error to give 


Et(kx.ky) 




-jkr 


2Ar 



^j(k^x + k^) 


dxdy . 


(2.1.9) 


Finally, a more compact representation, using vector notation for 
(2.1.9) is. 


--’o / *' '*5 • (2.1.10) 

s 

All the other equations could also be written in a similar form. This 
notation will be used in the following sections. 


2.2 The Antenna as a Receiver 


A general expression for the voltage at the terminal of an aperture 
antenna acting as a receiver is derived in this section. The antenna is 
represented by the physical aperture, S. The feed system is assumed to 
be accessible at some point so that voltage and current can be measured. 
We further stipulate that the wave guide or transmission line will sup- 
port only one propagating mode. Higher order modes are usually excited 
by reflection at the end of the feed section; but we assume that these 
are strongly attenuated before the point of measurement. 
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Figure 2.2. The antenna feed system. 


and are the voltage and current, respectively, measured when 
the antenna is used as a transmitter. and are similarly measured 
during reception. Directions of these quantities are defined in Figure 

2 . 2 . 

Let an incident electric field be denoted by E^. Its angle of 
arrival with respect to the antenna is measured in spherical coordinates 
as (0,<j)). The vector notation indicates the direction of polarization. 

The Rayleigh-Carson Reaction Integral theorem of mixed form [2,3] 
states that, 

''t'r - ''r't = f (2-2-D 
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and if 


Z, = Antenna Impedance 

a 


= Load Impedance 


Vt = It^a 


Vr = -Ir^l • 


Replacing fJg, Mg) by the equivalent current density 



( 2 . 2 . 2 ) 


Thus, the effect of the current distribution due to the antenna feed is 
to weight the incoming field and is included in (2.2.2). The actual 
implementation of the feed system, be it a horn at the focus of a para- 
bolic antenna or a system of wave guides connected to various points of 
the radiating surface, is immaterial to Equation (2.2.2). Hence, all 
classes of aperture antennas are included. 


2.3 Case 1: Plane Wave Incidence 


The general expression- for the antenna terminal voltage derived in 
Section 2.2 is applied to the plane wave incidence case. The identity 
of the transmitting far field and the receiving patterns are established. 

Consider a plane wave travelling toward the origin. Its angle of 
arrival is measured as shown in Figure 2.1. It may then be represented 
as 
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(2.3.1a) 



+ k^,y + k^z) 



-jkT 

e 


(2.3.1b) 


~ J k * ) * 

= p Aq e (2.3.1c) 

where p is the unit polarization vector of the plane wave. Aq is the 

constant amplitude of the wave. 

} 

Now 7^(^) was defined to lie only in the x-y plane. Therefore, 
from (2.2.2) 


V = 

/ (P) . E^(P) dS. 

(2.3.2) 

where 

It : 



1.-0 

(2.3.3) 

-j(k X + k^) 

(2.3.4) 

-jT.p 


= Aq e 


(2.3.5) 


is the field distribution of the incident wave over the aperture plane. 
Further, using (2.1.1) and (2.3.2) 
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V 


r 



J f(P) • ErC) 


Let the constant 


K 


V 





Then, 


V 


r 



dS 


(2.3.6) 


(2.3.7) 


or 

/ “jR” *P 
f(P) e " dS 


i .e. , from (2.1.6) 


(2.3.8) 


V^{ks) = K/qP . g(-¥^) . (2.3.9) 

Note that the sign of k in 7^ is negative for a wave travelling toward 
the origin of coordinates. Equations (2.3.9) and (2.1.8) are seen to 
have the same form. Therefore, the receiving pattern for plane wave 
incidence is identical to the far field transmitting pattern. 

Finally, using the approximations discussed in Section 2.1, the 
integral in Equation (2.3.8) may be carried out over the physical aper- 
ture only, so that 
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(2.3.10) 


/ - jk *P 

7(p) e * dS . 

S 

V^Ckg) gives the voltage at the antenna terminals due to a plane wave 
incident in a direction such that the projection of the propagation 
vector on the antenna aperture plane is 1<g. 


2.4 Case 2: Perturbed Wave Incidence 


The incident wave is now perturbed deterministical ly in amplitude 
and phase; and the antenna terminal voltage is derived. This is an 
intermediate step toward the randomly varying incidence case. 

Consider a wave of the form 
-jT<*r 

= A(r) e (2.4.1) 

incident on the aperture antenna. A(r) is the complex amplitude which 
is not constant over a plane transverse to the direction of propagation. 
This could be visualized as the result of an initially plane wave, with 
constant amplitude and phase in a plane transverse to its direction of 
propagation, which has undergone perturbations along its front. An 
example is the field near a small fixed scatterer. Using (2.3.3) - 
(2.4.1), 

-jl< •P' 

Hr) I = A(P‘) e ^ (2.4.2) 

I z=0 




= K, 


/ 


f(P) 


-jk-.p 
A(P) e 5 


dS 


(2.4.3) 


16 




PLANE WAVE ( PERTURBED WAVE 


SPACE I RANDOM 

' ATMOSPHERE 


J 

ANTENNA 


Figure 2.3. Perturbation of a propagating wave by the turbulent atmos- 
phere. 


2.5 Case 3: Randomly Varying Perturbed Wave Incidence 

The expression for the antenna terminal voltage will now be ex- 
tended to the case of a randomly varying incident wave. A plane wave, 
travelling from a spacecraft and incident on the earth's atmosphere, 
encounters inhomogeneities in pressure, temperature, and humidity (Fig- 
ure 2.3). Their influence is strongest in the lower troposphere due to 
changes in the refractive index. The effect of bulk refraction is to 
bend the propagating ray as a whole. This is a relatively slow pheno- 
menon and could be considered to be quasi-static over times on the order 
of several cycles of carrier frequency. It will not be considered in 
the following analysis. 

Inhomogeneities of smaller scale sizes are also present in the 
atmosphere. The wave, moving in a direction ¥ with respect to the an- 
tenna's frame of reference is increasingly perturbed as it propagates 
through these infiomogeneities. Further, as these fluctuations are ran- 
dom in space and time varying, the perturbation of the wavefront will 
also be random in space and time. However, it is reasonable to assume 
that the average direction of propagation of the wave as a whole is 
essentially unchanged. It is, therefore, convenient to separate out the 
basic propagating direction and the random fluctuations. 
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Let 


Ep = A(r,t) e 


-jk*r 


where 

_ x(r,t) + j%r,t) 

A(r,t) = Aq e 

The tilde denotes random quantities. 'K is a function of position and 
time. is the amplitude of the wave that would be present if there 
were no perturbations, x is the log amplitude of the wave and the phase 
fluctuations are given by 


Or 

The statistics of x and ’i' are directly related to the statistics 
of the medium. This will be considered in detail in Chapter Three. The 
processes are assumed to be stationary and, hence, ergodicity may be 
invoked. Therefore, we may restate the above equations as functions of 
space only as 

= A(r) e (2.5.1) 

_ x(r) + j^r) 

A(r) = Aq e . (2.5.2) 

Let the value of on the x-y plane be 


E^{P) = 1^(7) 


2=0 


'V. - jk .p 

= U(7) e ^ 


where 


xi* x(x,y,0) + j>(x,y,0) 

U(F) = 7^0 e 


(2.5.3) 
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(2.5.4) 



^ X(P) + j«(F) 
= 7^0 e 


P is the surface vector (x,y) as before. Evaluation at z=0 is under- 
stood and will be dropped from the notation in the following. 


The voltage at the antenna terminal when receiving such a wave 
would be 

^ J'fiP) •%{?) dS . (2.5.5) 

oo 

Once again, using the engineering approximations discussed in Section 

2 . 1 , 



S 


Equations (2.5.5) and (2.5.6) are the statistical extensions of (2.3.7) 
and (2.3.10). 


2.6 Fields in the Focal Plane of a Focusing System 


In the particular case of a focusing system, when the incident 
field is paraxial, the field in the focal plane of the antenna may be 
obtained fairly easily ([11], pp. 231-236). The results can also be 
directly applied to the optical telescope, which can be considered to be 
a focusing antenna with a uniform polarization insensitive aperture 
illumination function. Let 


P 

F' 

L((x,y) 


focal length of the antenna 

(x,y) be a point on the antenna aperture 

(x',y') be a point on the focal plane 

the magnitude of the incident field on the aper- 
ture plane 
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the magnitude of the field in the focal plane. 


Then, neglecting polarization, 

r ^ jf-(xx' + yy') 

U'(x',y') = / U(x,y) e ^ dS, (2.6.1) 

s 


Equation (2.6.1) is seen to be similar in form to the expression for the 
terminal voltage of a receiving antenna. Therefore, it serves to link 
the results obtained for an antenna to the optical telescope. 

Comparing (2.6.1) and (2.5.6), it is seen that if the following 
substitutions. 


k X'/f^ 

for 

k sin0 

cos4> 

k y'/fj 

for 

k sine 

sin<t> 

% 




U(x',y') 

for 

Vr 



are made, the results derived for the received voltage may be also ap- 
plied to the focal field. 

% 

Clearly, the description of the received voltage, V , or the field 
in the focal plane, U', in terms of statistical averages, requires know- 
ledge of the statistics of U(P). This is considered in Chapter Three. 

In summary, the theory of the receiving aperture antenna was con- 
sidered. The basic formulation was extended to describe the behavior of 
the antenna when a randomly varying perturbed wave is incident upon it. 
Finally, the field in the focal plane of the antenna was shown to be 
useful in relating the results obtained in this study to the optical 
telescope. 
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CHAPTER III 


PROPAGATION IN TURBULENCE 


3.1 Introduction 


The propagation of electromagnetic waves through the turbulent 
atmosphere has been the subject of intensive investigation over several 
decades, A vast literature exists, detailing several approaches to the 
problem. Some basic references [4-22], and works referred to in their 
bibl iography, may be consulted for comprehensive study. 

Chernov [4] and V.I. Tatarski [5,6] brought together much of the 
earlier work and also provided a broad theoretical basis for plane wave 
propagation through the turbulent atmosphere. This was extended by 
Schmeltzer [7], considering laser beams, to include spherical wave, fin- 
ite aperture cases. The Rytov method of evaluating the wave equation 
is used in these works. The range of applicability of the Rytov method, 
which requires, inter alia, that the magnitude of the fluctuations of 
the wave be 'small', sparked considerable debate [19]. Various mathemat- 
ical approaches, such as the mutual coherence theory [20], have been 
used to study the strong fluctuation case. 

In 1941, A.N. Kolmogorov proposed a theory of turbulence which is 
applicable to the earth's atmosphere [21]. This is now a generally 
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accepted description. Lee and Harp [22], used a more heuristic model of 
the atmosphere, i.e., a set of phase gratings, to obtain solutions to a 
wide range of problems. Ishimaru, in his books [8,9], provides a useful 
collection of these ideas, as well as some of his own, in a readable 
text. 


This chapter considers some aspects of the propagation of a plane 
wave in a turbulent atmosphere. A review of basic concepts is followed 
by comments on the Gaussian and Kolmogorov models for clear air propaga- 
tion. The simpler Gaussian model is also considered because it is more 
tractable, mathematically. The two models are then compared to obtain an 
understanding of the uses and possible errors associated with the Gaus- 
sian approximation. Some new expressions are also derived for log ampli- 
tude and phase, variance and correlation functions. 


3.2 Preliminaries 


The propagation of electromagnetic waves along terrestrial or 
earth-space paths in the microwave through optical portions of the spec- 
trum is most influenced by that portion of the earth's atmosphere called 
the troposphere. This denotes that part of the atmosphere extending 
from the earth's surface up as high as 17 km. The refractive index, n, 
of the atmosphere is approximately one and is given by 


n = 1 + An 


(3.2.1) 


and the refractivity, N, is given by 


N = (n-l)xlO® = 

T 


+ 4810 e^ 


in the microwave region and 


(3.2.2) 
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N = (n-l)xlO® = ZZjlLZ 

T 


(3.2.3) 


in the optical region. Here, N is on the order of 300 and 

P = Barometric pressure in millibars (1 mm. Hg = 1.3332 millibars) 
e = Partial pressure of water vapor in millibars 
T = Absolute temperature, degrees Kelvin. 

The typical decrease of n with increasing height causes any pro- 
pagating ray to bend downwards. This can be eliminated in calculations, 
under conditions of standard refraction, by considering the earth to 
be a sphere of radius equal to 4/3 the physical radius, or, approximate- 
ly 8479 km. This approximate correction will be used throughout this 
work. 


The pressure, temperature, and humidity also vary randomly from 
place to place and with time at any place due to atmospheric turbulence. 
The resulting random variation in the refractive index perturbs the 
propagating electromagnetic wave causing amplitude scintillation and 
phase fluctuation. Other atmospheric effects such as absorption and 
depolarization have been extensively treated in the literature and will 
not be considered here. 

The correlation function, B^(Fpr 2 ), of a multi-dimensional random 
field, f(r), is defined [6] as 


BfCrpr^) = <[f(rj) - <f(rj)>] [f(r 2 )*- <f(r 2 )*>]> . (3.2.4) 

The random process is said to be homogeneous if it has a constant 
mean and its correlation function is unaffected by simultaneous transla- 
tion of r^ and r 2 by the same amount in the same direction. It is the 
generalized concept of a stationary process. The correlation function 
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of a homogeneous random field depends only on the difference r = Tj-r 2 , 
so that 

^f('"l»'"2^ “ B^(rj^-r2) = B^(r) . (3.2.5) 

A homogeneous random field is said to be isotropic if B^(F) depends 
only on r = [Fj, i.e., only on the distance between the observation 
points. 


Bf(ri,r 2 ) = B^(r) . (3.2.6) 

Many meterological phenomena are not properly described by homo- 
geneous, isotropic, random processes. However, in many cases a new 
process f(F^) - f(F 2 ) can be formed which is well behaved (homogeneous). 
f(F) is then called a locally homogeneous process. Its correlation 
function is no longer a function of F 2 -F 2 alone; but also depends on 
Fj and F 2 , individually. Therefore, the structure function 

Df(F) = <lf(F^+F) - f(F^)|2> (3.2.7) 

is a more useful parameter of such a field. It is a function of F 2 -F 2 
only. Furthermore, if D^(F) is a function of r only, the process is 
called a locally homogeneous isotropic random field. 

It should be noted that the structure function exists for homogen- 
eous as well as locally homogeneous random processes. A unique correla- 
tion function, on the other hand, exists for homogeneous processes only. 

The three-dimensional correlation function, B^(F), and its spectrum 
•lyCic) are related by the Fourier transform [23], 

1 r - 

1'f(<) = — ^ / B.(r) e dF (3.2.8) 

^ (2tt)3 J ^ 
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(3.2.9) 



is called the spatial frequ^epcy, or the turbulence wave number. £ is 
the size of turbulence eddies. 

Also, for homogeneous isotropic turbulence, 

D^(r) = 2 [8f(0) - Bf(r)] . (3.2.14) 

Since B^(0) = if b^(r) is the correlation coefficient, 

D^(r) = 2o?(l - b^(r)) . (3.2.15) 
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Finally, for the correlation function to exist in general, the 
spectrum, $^(k), must be regular. However, the structure function can 
exist even if the spectrum has a singularity at the origin. This is a 
particularly useful property of the structure function. 


3.3 The Gaussian Model of Atmospheric Turbulence 


In this model, the turbulence is assumed to be homogeneous and 
isotropic. Consequently, the correlation function of the refractive 
index fluctuation exists and is a function of spatial separation, r, 
only. It is defined by 


B 


ng 


(r) = a 


2 

n 



(3.3.1) 


where 

p 

at = variance of n 
n 

= correlation length . 


The correlation coefficient is then 

.ri 

in 

b„g(r) = e " . (3.3.2) 


Thus, Z is the separation at which the correlation coefficient falls to 
1 ^ 

- . Hence, for the Gaussian Model, from (3.2.11) through (3.2.13), 





(3.3.3) 
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(3.3.4) 


= 2a^,(l - b„g(r)) 


(3.3.5) 


Note that (3.3.3) may also be written as 







where 


(3.3.6) 

I 



This representation is useful when comparing the Gaussian spectrum with 
the Kolmogorov and Von Karman spectra. 


3.4 Modified Kolmogorov Models 

The turbulence spectrum of the Kolmogorov model is divided into 
three regions by two scale sizes: the outer scale of turbulence, Lq, and 
the inner scale (or microscale) of turbulence, Hq [25,26]. Let Jibe the 
size of turbulence eddies. 


1. Input range Lg < £ , ^ ^ 

Large scale atmospheric characteristics such as wind form turbulence 
eddies of this dimension. Thus, energy is input to the turbulence system 
from the thermal and kinetic energy of the atmosphere. The process is. 
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in general, anisotropic and varying, depending on climatic conditions. 
Consequently, it has not yet been described successfully in closed form. 
The spectrum is, therefore, undefined in this region. 

The value of Lq may vary widely. Close to the ground, it is usual- 
ly taken to be equal to the height from the ground. In the free atmos- 
phere it may be in the range of 10-100 m or more, differing significantly 
in different regions of space and depending on the turbulence inducing 
mechanism. 


O 0 

2. Inertial range ^0 ^ ^ *"0 ^ ^ ^ 

The eddies formed in the input range are unstable and fragment 
into smaller ones. These break up as well, continuing in this manner 
and causing energy to be distributed from the small to large turbulence 
wave numbers. There is very little energy loss in this process. Kolmo- 
gorov has shown that the spectrum in this region follows the law, 

_11 

^nk^*^^ K (3.4.1) 

and may be considered to be isotropic. Most cases of microwave propaga- 
tion are affected predominantly by this region of the wave number spec- 
trum. 


3. Dissipation range A < Jl«, T* ^ 

U X.Q 

The energy in the turbulence, which was transferred through the 
inertial subrange, is dissipated through viscous friction by very small 
eddies. An accurate characterization of the spectrum is not known. It 
tends to zero rapidly and is frequently assumed to be of the form. 
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2 

s_ 

2 



(3.4.2a) 


with 


'^m 


5.92 


(3.4.2b) 


2.Q is variously assumed to lie in the range of 1-10 mm. The composite 
spectrum, in the region other than the input range has been expressed 
by Tatarski as [6J 




0.033C^ K 


11 

.3 


'm 


(3.4.3) 


2 

where is the structure constant which will be considered in greater 
detail below. 

The spectrum of the atmospheric turbulence is locally homogeneous. 
Hence, (see Section 3.2), a correlation function does not exist. However, 
the structure function exists and is related to ^f,(<) by the relationship 
(3.2.8). 

In the inertial subrange, which is approximately isotropic, if we 
approximate the full spectrum by (3.4.1), we can write the structure 
function of the Kolmogorov spectrum, using (3.2.13), as 

2 

^nk^*^^ = Lq » r » . (3.4.4) 

This defines the structure constant in this case. 
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Von Karman Model 


In order to obtain a full range solution using the Kolmogorov 
spectrum, the input region needs to be defined. One approximation which 
eliminates the singularity at the origin in (3.4.3) and gives a homo- 
geneous isotropic spectrum is [27,28]. 

11 ■ 2 


where 



(3.4.6) 




(3.4.7) 


Also, since Hq is of the order of a few millimeters, is very large. 
Assuming that -> <», the corresponding structure function is, [28], 



(3.4.8) 


where Kj^(x) is the modified cylindrical Bessel function. (3.4.8) should 
3 

correspond to (3.4.4) for small r. Hence, 
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.2 

p2 % , Q, . 32 

Cp -V 1.91 Lq Op 


m 


(3.4.9) 


Comparing (3.2.15) and (3.4.8), we obtain the correlation function 


B (r) 
nv' 


2 2'" 
^n nv 


Wo. 


(3.4.10) 


Therefore, when so that (3.4.9) is applicable, 




£. ± 


r ' 


(3.4.11) 


and 




0.033Cp(K^ + <^) ® e 


(3.4.12) 


3.5 Averages oF the Components of the Received Field 

The random field in the aperture plane of a receiving antenna 
was described in (2.5.4) as. 


'V' x(p)+jX^) 


U(P) = Aq e 


(3.5.1) 


In this section we obtain expressions for the variances, and o^. of 

'V, 'V, X 4' 

the log amplitude and phase x and i};, respectively. 


Let, 
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L = the path length travelled by the wave in the atmosphere, 
p = the distance between two points in a plane transverse to 
the direction of propagation. 

Then [29], assuming a homogeneous isotropic medium, it can be shown that 
for any quantity, x, the transverse correlation function is 

00 

B^(p) = 2TT^k^L J fj^(ic)$^(ic)jQ(Kp)< d< , , (3.5.2) 

0 

where 

fj^(ic) = the spectral filter function of x. 

h 

This is so named since it, together with the zero^" order cylindrical 
Bessel function, filters the spectrum, to give the correlation 

function B^. It is obtained from (3.5.1) using Rytov's first approxima- 
tion. Then, for a plane wave, the log amplitude and phase filter func- 
tions are [29], 


fx(K) 



(3.5.3) 





(3.5.4) 


The transverse correlation functions may be defined as. 


By(p) = Oyb (p) = the log amplitude correlation function, 

A A A 

B^(p) = a^b^(p) = the phase correlation function. 
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o 

c and b denote the variance and correlation coefficients in the trans 
verse plane, respectively. ' The following results are derived in Appen 
dix A2 and simply stated below. 


Gaussian Spectrum 


From Appendix (A2.1), 




(3.5.5) 

(3.5.6) 


where the Gaussian wave parameter is 


. (3.5.7) 


Also, if we define the wave correlation function as 


B^(p) = Bj^(p) + B^(p) , 


then 



all„A e 



(3.5.8) 


(3.5.9) 


33 







■^g ] 


9 



Wote 






+ o^bJP^ 


'fV 




O-S'iol 

(3.5-"^^^ 

(3. 5.^^'* 


(,3.^ 


.I3a) 


(3.^ 


.13'3') 




, . .2a - “ “ 

.2 a - 


= =2.a ' ■ 


.476 


5ei of «>'®® 

,3 2 IS') ’ 

fro^a ,4Ann 


in 


(3.5 




to 5e 


(3.5*^^''^ 

one-V^a''^ 


tV\e 


corresp 


ond^^9 


(3.5 






0 


34 



Von Karman Spectrum 


From Appendix A2.2, 


BJP)' 

^(^>1 


5- 

0.033A^C^k^L 


Y. -fiHjCrf) 


2p 


p=0 <P') 


✓ 

r( P+1 )'i'( P+1, p+^,Zq) 

\ 






(3.5.14) 

(3.5.15) 


where 


Z = Zq - jWo 


Zn = 


1 


<2,2 

m'-O 


= 


kU 


"m 7^ k 


(3.5.16a) 


(3.5.16b) 


Wq and Wjjj may be called wave parameters for the Von Karman model. Next, 


5 

B^(P) = 0.033it2 2Lhlkh 


^ (p!)^ 

p=0 



r(p+l)y(p+l. p4>l, Zq) 


(3.5.17) 
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' 5 


■ .11 ■ 


’ = 0.033tt 3 I L^C^k^L • ' 


(-1) - J“o) 


/ 

0 

L \ /J 

> 


(3.5.18) 

(3.5.19) 


5 

2 „ 12 , 3r2, 2 , 

% = 0-033’' -3 LgC^k L 


(3.5.20) 


Kolmogorov Spectrum 

This spectrum is not homogeneous and isotropic. Therefore, the 

correlation functions are not defined. The singularity at the origin of 

the extended spectrum (3.4.3) results in being singular. 

However, f y( '^) '*’ni,( '^) is finite for all k. Thus, only may be calcu- 
X ni( X 

lated. This is obtained in Appendix (A2.3), 


= 0.033if 
X 


5 


2 6 rn 


11 . ^11 

ifci + w3)12 


tan'^W ) 
m/ 


m 


(3.5.21) 


Note that if W, 


m 


7 n 

1, this reduces to the common k L form. 
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3.6 Some Properties of Propagation in Turbulence 

Some general properties of propagation in homogeneous turbulence 
are discussed in this section using the concept of filter functions. A 
detailed treatment may be found in [9] and is not attempted here. 


Consider a medium having a Gaussian spectrum as an example. Fig- 
ure 3.1 shows such a spectrum together with examples of the amplitude 
and phase filter functions, ^nd f^(<). 

Equation (3.5.2) shows that the integrand of the expressions for 
amplitude and phase variance include the product of the turbulence 
spectrum and the corresponding spectral filter function. Thus, it will 
be seen from Figure 3.1 that the phase fluctuations are affected by all 
eddy sizes and especially by those in the input region. On the other 
hand, for the log amplitude fluctuations, setting 


yo = 1 

to determine the cut-off point yields 

2 , 

sin(^] 0 

so that the smallest < for which f {<) = 1 is 

X 



For smaller wave numbers, the amplitude filter function tends 
rapidly to zero, while for larger wave numbers, i.e., smaller eddies, 
the value of the spectrum is smaller. Hence, the major contribution to 
log amplitude fluctuations comes from eddies having sizes of about 



i.e., comparable to the size of the Fresnel zone. 
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Figure 3.1. Atmospheric spectrum and filter functions. 
= 46 m, F = 30 GHz, L = 100 km. 
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Case 1: Geometrical Optics Region (Wg << 1) 


W 


9 



Then 


This case is illustrated in Figure 3.2. 
From (3.5.10) and (3. 5, .11), using 


tan-^Wg S: Wg - 

3 


^ 0^11 
2 n n"^ ^ 


1 + 1 - 




SO 


X 3 n J^3 
n 


(3.6.1) 


0^ - •./tT" a^z k^L 
^ V n n"^ ' 


(3.6.2) 


2 . . q 

Thus, 0 ^ is independent of frequency and is proportional to L"^. 


2 '' ? 
However, is proportional to k and L. 
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Figure 





3.2. Spectral filter functions for the geometrical optics region. 
= 100 m, F = 90 GHz, L = 0.5 km. 
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Case 2: Long Path Region (W 1 ) 


In this case 




This is illustrated in Figure 3.3. Again, from (3.5.10) and (3.5.11) 
and using 


tan" W. 


-9. << 




(3.6.3) 


Now both 0 ^ and are approximately equal and are proportional to 
and L. 

2 2 

remains proportional to k L in both the geometrical optics and 
long path regions. 


Though the Gaussian spectrum was used here as an example, it is 
important to note that the same behavior is true for any homogeneous 
isotropic spectrum. 


In a medium having a Von Karman spectrum, 


w << 1 
”m ^ 


denotes the geometrical optics region, with 


-V jj2 .V, |^2|^ 


1 << Wr 
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Figure 3.3. Spectral filter functions for the long path region. 
= 10 m, F = 0.1 GHz. L = 300 km. 
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represents the long path region, where 
'V' k^L 

Also, 


represents the intermediate region in which the propagation characteris- 
tics are dominated by the inertial subrange. Here 

'V k\ ^ . o| -V- k^L . ! 

If the medium is represented by the extended Kolmogorov spectrum 
alone, there will be no long path region. Then, the wave parameter, W^, 
will divide the spectrum into the inertial subrange and the geometrical 
optics region only. The remarks on these regions, made above, would 
still apply. 


3.7 Comparison of Spectra 

Some main features of the Gaussian and Von Karman spectra are 
compared in this section and are related to each other. 

Figure 3.4 shows typical correlation functions for both the Gaus- 
sian and Von Karman spectra, (3.3.1) and (3.4.11), respectively. The 
main parameters in this plot are 

= 40 m 

Lq = 52 m 
= 10 mm. 
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Figure 3.4. A comparison of correlation coefficients of the Gaussian 
and Von Karman models, 
tf, = 40 m, Lq = 52 m, £q = 10 mm. 
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In an attempt to find a meaningful comparison, the parameters of 
the Von Karman spectrum were chosen so that both spectra showed the same 
correlation length, defined as that length at which the correlation 
coefficient drops to 1/e. The correlation length of the Gaussian spec- 
trum is, of course, It was found, somewhat surprisingly, that over 

a fairly wide range of values of the simple relationship 


= 1.2 to 1.4 (3.7.1) 

n 

appears to hold if the correlation lengths are forced to be equal. Under 
these conditions, the Gaussian correlation function shows a higher value 
than the Von Karman spectrum over small separations. However, for fair- 
ly wide separations, the Gaussian function cuts off sharply while the 
Von Karman spectrum exhibits a higher, although insignificant, correla- 
tion. 


Figure 3.5 shows the corresponding spectra, given by (3.3.3) and 
(3.4.12), respectively. The most clearly defined difference is the much 
higher value of the Von Karman spectrum throughout the input region. 
There is a short region where the spectra cross over. The slopes are 
fairly close in a portion of this region. For high values of k, the 
Von Karman spectrum, once again, is higher. In light of the preceding 
discussions on amplitude and phase filter functions, these characteris- 
tics imply that the variance of a propagating wave will be higher in the 
long path and the geometrical optics regions of the Von Karman spectrum. 

Another comparison was also made. The expressions for the wave 

variances of the two spectra are given by (3.5.12) and (3.5.20). Since 

2 

both quantities show a k L dependence, the condition under which they 
are identical is easily found. If 

a = a 
wg wv 
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Figure 3.5. A comparison of spectra of the Gaussian and Von Karman 
models. 

= 40 m, Lq = 52 m, = 10 mm. 


/ 
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-f <•1^/'- - 0-033"^c2(l|)L3k2L . 

Then, using (3.4.9) we obtain 

Lq = 1.2 . (3.7.2) 

Comparing (3.7.1) and (3.7.2), it is clear that (3.7.2) is a reasonable 
basis for matching the Gaussian and the Von Karman spectra. 

The discussion so far has not included the variance of the refrac- 
tive index fluctuations. This may be accomplished by using (3.7.2) 
together with (3.4.9) to define an equivalent structure constant, C^g, 
for the Gaussian spectrum. 


2 

C^g = 1.91 (1.2 (3.7.3) 

This quantity, by itself, is meaningless. It does serve a useful pur- 
pose, however, as an aid for comparing the two spectra. This is especi- 

2 

ally useful as most experimental and theoretical works describe C and 

2 " 
not o„. 

n ’ 

For the purposes of engineering calculations, the Gaussian spec- 
trum provides a reasonably good approximation to the Von Karman spectrum. 
Its advantages for analytical purposes are many. Its major sources of 
possible error, if borne in mind when interpreting the results, are 
tolerable. Consequently, the analysis that follows will assume a homo- 
geneous, isotropic Gaussian spectrum, given by (3.3.1) and (3.3.3). 
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3.8 Moments of the Aperture Field 


The moments of the aperture field of an antenna in a turbulent 
medium are given in this section. The following assumptions are made: 

1) The turbulence spectrum is assumed to be homogeneous, isotropic, 
and Gaussian. Its correlation function is, therefore, also Gaus- 
sian. 

2) The average direction of propagation of the wave, as a whole, is 
unchanged. Also, there is no depolarization of the wave. 

3) The antenna is relatively large in terms of wavelength. Hence, 
the beamwidth is fairly small, i.e., on the order of a few degrees 
at most. The side lobe levels are also assumed to be low. 

The third approximation above implies that if the analysis is 
confined to the main lobe, even down to the first null, the aperture 
plane is nearly transverse to the direction of propagation. Thus, cor- 
relation functions and variances in the aperture plane are substantially 
the same as those derived for the transverse plane. This is referred to 
as the paraxial approximation. 

Under these assumptions, if U(P) is the random field in the aper- 
ture plane, then from Appendix A3, letting 


^0 " 




Ox + jOxO 




(3.8.1) 


we find that 


< U(P) 


> 



2 w 


and 


(3.8.2) 
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< 


(3.8.3) 


U(l"l) U{J^) > = |U/ e % 


oh (p) 

w n'^ ^ 


where 


p = 



and 


Uq = pUq . (3.8.4) 

p is the polarization vector as before. 

It is interesting, as pointed out by Shifrin [11], that under the 

% *\t if 

assumptions made for this analysis, <U(P 2 ) U(P 2 ) > remains the same 
whether the amplitude and phase perturbations are mutually dependent or 
not. Therefore, the mean power pattern and other parameters associated 
with the power pattern are also independent of the nature of the interac- 
tion between the random amplitude and phase, if any. 


3.9 Incident Power Density 

The expected value of the time average power density of the inci- 
dent field in the antenna aperture is 


S. 

me 


> - 


1 

7 


Re 


< i(F) 


i(F)*> 


(3.9.1) 


where H^(P) is the random incident magnetic field in the aperture plane. 
If the perturbations are such that the incident wave can be considered 
to be locally plane, then taking Hq as the characteristic impedance of 
free space. 
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< 


'V/ 


> 


Ejp) • r (p)* 

< c . > = 5 : 

me 2 hq 


(3.9.2) 


Therefore, from (3.8.3), 


< S • > 

me 




(3.9.3) 


Fried [34] has shown from energy eonsiderations that, for a plane 
wave propagating through the turbulent atmosphere, assuming no out of 
path seattering. 





(3.9.4) 


This implies that the power density, < is eonstant. These as- 

sumptions will not, however, be made in this work. Uq will be left in 
its full form. 

Summarizing, the main thrust of this ehapter was to find moments 
of the field in the aperture plane of the antenna. To obtain these, the 
general characteristics of some models of atmospheric refractive index 
fluctuations were studied. The averages of the components of the re- 
ceived field were also obtained for these models. The Gaussian model 
for homogeneous turbulence was related to the Von Karman model and was 
found to be a reasonable approximation for engineering calculations. 
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CHAPTER IV 


THE RECEIVING SYSTEM 


4 . 1 General Comments 


The results derived so far concern the effect of the atmosphere 
and the antenna only. In this chapter we examine some effects of the 
receiver on the reception of the scintillating signal. 

The receiver will be defined as that section of the receiving 
system which follows the antenna. Its role is to extract information 
from the incoming signal. To do this it may use, in some manner, the 
amplitude, phase, frequency, or a combination of these parameters, of 
the antenna output. The type of processing has a profound influence, 
not only on the statistics of the final output of the receiver, but 
also on the perceived statistics of the incoming signal as measured 
through the receiving system. 

In this study, the final output of the receiver will be called 
the output signal, while the incoming signal as measured through the 
receiver will be termed the measured signal. It is important to realize 
that in general the output signal, the measured signal, and the incoming 
signal, as defined above, will not be the same. The measured signal is 
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deduced from the output signal using the inverse of the receiver charac- 
teristics. However, due to the processing done by the receiver, some 
of the attributes of the incoming signal may be irretrievably lost. 
Therefore, the measured signal may be only a partial description of the 
incoming signal. An obvious example is that of the frequency modulation 
(FM) receiver, which eliminates amplitude information from the incoming 
signal through the use of limiters. 

This study will concentrate on receivers which do not destroy 
amplitude information. Expressions are derived for the received signal 
voltage level, non-fluctuating power or dc power, total power, normal- 
ized variance, and the degradation of the dc power, as measured by the 
receiver. For this purpose, the antenna will be assumed to be conjugate 
matched to the receiver in the rest of this study. 

Further, the average effective aperture, gain, and gain degrada- 
tion of the antenna are also obtained. These parameters of the antenna 
do not depend on the receiver used. 

The receiver will be divided into one of two broad classes as 
follows: 

a) Synchronous receivers. 

b) Asynchronous receivers. 


Synchronous Receivers 

A receiver using phase locked loops in its demodulating chain 
falls in the first category, as does any receiver using the phase inform- 
ation contained in the signal. It should be noted that an amplitude 
modulation (AM) receiver using a synchronous detector belongs to this 
group as well. Such a synchronous receiver's output amplitude, is 
affected by the amplitude scintillation and phase fluctuations of the 
incoming wave. Therefore, letting the tilde denote random quantities as 
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before, we may write 

'V 'V/ 

Vpg = function of (V^) . 

A receiver that employs limiters on the signal amplitude and uses 
the phase information only, is affected by phase fluctuations. If 
is the measured phase at the output, 

(4>Qg) = function of (Arg (V^)) . 

Consider a synchronous amplitude modulation receiver. Let be 
the received signal voltage as measured by this receiver. Then, since 
no amplitude or phase information is destroyed by the processing circuit- 
ry. 


V = RX' (V ) = V 
syn s ^ os' r 


(4.1.1) 


RX“^ is the inverse of the synchronous receiver transfer function and 
represents the receiver calibration curve. 

Next, consider the average value of 


< V 

syn 


> 


% 




(4.1.2) 


From (2.5.6), the mean received voltage is 


^ V 

syn 


> 



/ 


< U{P) > 


?(P) 


dS 
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Figure 4.1. The received signal measured by a synchronous receiver. 

n’ 


'\j 

< V 


syn 


> - 


Vp > = Ky J (p.Pf) < U(P) > f(P) e " dS . 


(4.1.3) 


The mean squared received voltage is 




JJ < u(F^)u(P2)* > 

s s 


f(Pl)f(P2) e s i ^ ^ 


We define the normalized variance as 
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(4.1.5) 


Asynchronous Receivers 


A receiver using a square-law detector to process amplitude inform- 
ation belongs to this category. This type of receiver does not use the 
phase of the incoming signal. Hence, it is sensitive to amplitude scin- 
tillation of the antenna voltage only. Any fluctuation of the overall 
phase front of the incoming wave will not be detected by the receiver, 
except in the sense that it affects the amplitude of the antenna output. 
If Vqq is the receiver output voltage, we may write 

= function of (Iv^l) . 

In the particular case of a receiver using square-law detection. 


Consequently, fluctuations of as measured by this receiver will be 
less than that measured by a synchronous receiver. This is an important 
difference. 

i 

Let be the received signal as measured by an asynchronous 

amplitude modulation receiver. Then, 


'^asyn *^^a ^'^oa^ ~ ^'^r^ 


(4.1.6) 


where RX~^ represents the calibration curve of the asynchronous receiver. 
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Figure 4.2. The received signal measured by an asynchronous receiver. 


Considering average levels. 


< 

*asyn 


> = < 



yy* < u(Pi)u(P2)*> ip-Pfi^ 

s s 


* -Jks*(Pl-P2) 


as before. 


It is not possible 
formation in the output. 


% 

to obtain < since there is no phase in- 

However, < Iv^l ^ may be obtained as follows 


'Xj 

^ V 

asyn 


> 



(4.1.7) 
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The normalized variance of the signal, as measured by an asynchronous 
detector, is 


2 < lv;|2 > - < Iv.l >2 

^'asyn " ' < i >2 

r 


(4.1.8) 


4.2 Average Received Voltage (synchronous receiver) 
From (4.1.2) 

'V 'V _ 

<v > = < y (k ) > 

syn 

and from (3.8.2) and (4.1.3) 

An^ 

% 2 w 

^ V ^ - Uq e Kv • y f(P) e dS 

S 


Let 



-jk -P 

f(P) e dS 


and 


(4.2.1) 


(4.2.2) 


(4.2.3) 


max ~ j^V^^s^jmax (4.2.4) 

over all 

Define 
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(4.2.5) 


9v(ks) 


V^s) 


I 


V max 


< 1 


and 


''dO = Wy max 


(4.2.6) 


Then, 


< V 

syn 


= V 


dO 


'2 w 
e 9v 


(k^) 


P'Pf 


(4.2.7) 


Vdo is the measured signal voltage when the wave is incident along 
the main beam axis, with no polarization mismatch. 9 y(l< 5 ) is the normal- 
ized plane wave voltage (or electric field) pattern function. Therefore, 
Equation (4.2.7) shows that though there is a reduction in the received 
signal voltage, there is no change in the field pattern function when 
the measurements are made using a synchronous AM receiver. 


4.3 Total Received Power 


The total received power, assuming a matched load is 



< 


lv^(ks)l^ > 

2|Z/ 


(4.3.1) 


is the load resistance, being the real part of Z 


1 


= Z 


* 

a' 


From (3.8.3) and (4.1.4), 
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< 


(4.3.2) 


IV^(ks)!^ ^ = Ip-Pfl^ |Uq 1^ e iKyl^ Ip(l<^) , 


where 


Ip(ks) - JJ 


%bn(Pl.P2) 


S S 


_ * -3ks*(Pi-P2) 
f(l^l)f(P2) e dS^dS^ . 

(4.3.3) 


Consider the unperturbed plane wave incidence case, i.e., = 0. The 

integral becomes 


Vo'// 




dSj^d$2 


S S 


f(Pl) e dS^ 


/ f(P2) e 


... 


dS, 


I 


-jkc- P 

f(P) e dS 




from (4.2.3). This is the plane wave power pattern. Further, let 


I 


PO max 




2 

max 


(4.3.4) 
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and 


9oC<s) = 


^PO max 


9v(ks)|2 


(4.3.5) 


go(ks) is the normalized power pattern function of the antenna when 
receiving an unperturbed plane wave. Then, (4.3.2) becomes 


^IVy.(k5)l^>= lp*p^l 


IUqI^ 


^^Vmax* 


9o(ks) e 


-°w 


PO 


(k^) 


Now define 


9d = e 



Ip(ks) 

Ipo(ks) 


(4.3.6) 


Using (4.2.6) we obtain 


< iV^df;)!^ > = |V^g|2 90(l<s)9d iP-Pfl^ • <'*•3-^) 


Therefore, the received total power 


'’r ' '’R09oC‘s)9d IP'Pfl 


where 


'dO' 


RO 


2|Z 


"2 ''1 


(4.3.8) 


(4.3.9) 
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t; 


Pro is the maximum received power for plane wave incidence with no polar- 
ization mismatch, g^j is called the gain degradation factor. It accounts 
for the reduction of the measured total power (and gain) when the inci- 
dent wave is perturbed by turbulence. 

A method of evaluating (4.3.6), which also helps to bring out its 
physical properties is by using the relationship 

^ (°^b )™ 
m=l 

in Equation (4.3.3). Then, 



IpCkj) = I 


PO 




z 

m=l 


j2m 

w 

m! 


'Pm 


(kj) 


(4.3.11) 


where 


•jks-(Pi-P2) 


^Pm^'^s^ " // f(^l)f(P2)* e ' ’ " dS^dS^ . 


S S 


Let 




^PmO<s) 

Ipo(^s^ 


(4.3.12) 


(4.3.13) 


where will be called the correlation integral. It will be shown 
later that this is a function of the degree of correlation of the 
wave over the antenna aperture. 
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Thus, we may write (4.3.8) as 


-a 




1 + 


" 2m 

y ^ 

L. m! 


m! ^Cm^*^s^ 


m=l 


f Ip-Pfl^. (4.3.14) 


Now it is seen that 


„2 “ „2m 

-a — , a 


n W V" W T (T \ 

9d = ^ L ST 'erne's) • 


m=0 


(4.3.15) 


'cod's) “ 1 • 


(4.3.16) 


Finally, we let 


PROgoC<s) e 


w 


<v 

< V > 
syn 


2|Z^| 


2 *^1 


= P 


dcs 


(4.3.17) 


which we identify as the measured dc power, Pj^-g. using the synchronous 
receiver. It corresponds to the coherent power in the wave, as defined 
by Ishimaru. The second term in (4.3.14) must, therefore, represent the 
fluctuating power. Hence, 


Pn = P 


dcs 


+ P 


f 1 


(4.3.18) 


where 


'"fl ■ ^RO^O^^s^ 



m=l 


'cm(^s) 


(4.3.19) 


The fluctuating power, P^^ corresponds to the incoherent power in the 
wave as defined by Ishimaru. 
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4.4 Signal Variance (synchronous receiver) 

The normalized variance of the measured signal, from (4.1.5) and 
(4.3) is 


’’r * 


I I 

< V > 

' syn ' 


^syn^*^s^ 


2|Z 


1 ' 


'V/ I 

< V > 
syn ' 


2lZ 


1' 


(4.4.1) 


~o. 


9o(ks)gd - e Igv0<s)l^ 


W 


IgyCkg) 


Therefore, 

^syn^^s^ ® 9jj -1 • 

If (4.3.15) is used, we have 
» 2m 

^syn^^s^ " X! inf" ^Cm^^s^ • 

m=l 

Final ly, 

^syn ^°giO ^syn * 


(4.4.2) 


(4.4.3) 


(4.4.4) 
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4.5 Average Signal Level Degradation (synchronous receiver) 


The average signal level degradation may be defined as the reduc- 
tion in the measured dc power compared to that in the absence of 
turbulence. 


^^dcs ~ P 


dcs 

dcs 


(dB) 


No fluctuations 


(4.5.1) 


For the synchronous receiver, from (4.3.17) 
•"dcs " ''roSo^^s) ® 


dcs = -10 '“310 “ " 

(dB) 

(4.5.2) 

dcs = “w 

(dB) . 

(4.5.3) 


It should be noted that is independent of the antenna para- 

meters. This is because it is a measure of the reduction in the true 
coherent power in the incoming wave. This is clearly a function of the 
turbulence only. 
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4.6 Average Received Voltage (asynchronous receiver) 


% 
< V 




asyn 


(4.6.1) 


From Appendix A4, for small fluctuations, 

f 


^''asyn^"l''do‘ 'gv^^s)! 


4o 


gd]"^[^ ■ 23d ® ^]^|P'Pfl • 

(4.6.2) 


The steady state or dc power of the measured signal is then 

.2 


% 
< V 


asyn 




Therefore, 


’’dca ^ ^R03 q^'^s^ ‘ 


4o2 

. - 2g, * e *] 


< 

> 


|P*Pfl^ 


(4.6.3) 


This should be compared with the dc power for the synchronous receiver 
(4.3.8). Note that whereas the latter measures the true coherent power 
in. the wave, the asynchronous receiver measures the dc power in the 
envelope of the antenna output. It is insensitive to fluctuations in 
the phase of the incoming signal. 


4.7 Signal Variance (asynchronous receiver) 


s 


2 

asyn 




(4.7.1) 
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0 


s 


2 

asyn 




- 1 


4o 


+ e 


Also, 

^Lyn ' 1° '“9l0 ^Lyn • 


(4.7.2) 


(4.7.3) 


Equation (4.7.2) should be compared with (4.4.2) for the synchron- 

p 

ous receiver. s:^.,„ is a measure of the fluctuating power in the 
dsyn 

envelope of the antenna output signal. 


4.8 Average Signal Level Degradation (asynchronous receiver) 


Using (4.6.3) and noting once more that 


dca 


= '’r090<''s) 


No fluctuations 


AP 


dca 


-10 1 


O9l0 P 


dca 


deal 


I No fluctuations 


(4.8.1) 


The measured signal level degradation is 


4Pdca - -1° 


4=2,. 




> (dB) 


(4.8.2) 


APdea reflects the phenomenon that as the turbulence increases, a portion 
of the average power in the antenna output is converted to modulation 
side bands. 
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4.9 Average Effective Aperture 


The average effective aperture, Ag, of the receiving antenna is 
defined as 


fl = < T ot al Received Power > 

"e < Incident Power Density > 


(4.9.1) 


From (3.9.3) and (4.3.8), with matched polarization, and load. 


/ 2 


A = 


2|z,r 


2n, 


Therefore, 




(4.9.2) 


where 



I 


PO max 


(4.9.3) 



(4.9.4) 


AgQ is the maximum effective aperture of the antenna, when receiving an 
unperturbed plane wave. 
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4.10 Average Directive Gain 


Define the average (degraded) gain as 

Gj = • Average Effective Aperture . 

X 

From (4.9.2), we may write 

Sdli's' = ®0 max 9offs>9d • 

where 

p 

0 max 

- ll Ik |2 T 

- ^2 Vo max ' 

Furtherj if we expand g^j, as in (4.3.15), 

2 “ 

*^d^^s^ " % max 9o^^s^ ^ ^ Ys 

m=0 

From (4.10.4), we may define 


2m 

% 

m! 


Cm 


(ks) 



" ^0 max 90^'^s^ * 

Gq is the maximum directive gain for unperturbed 
incidence. C<Q(k^) is the plane wave incidence, power gain 
Hence, 


(4.10.1) 


(4.10.2) 


(4.10.3) 


(4.10.4) 


(4.10.5) 

plane wave 
function. 



4.11 Gain Degradation 


The gain degradation of the antenna in a direction is defined 
as 


AGCkg) = Go^'^s^dB " ^d^*^s^dB* 


(4.11.1) 


where 


®dB ~ ^®9l0 ^ * (4.11.2) 

Then, from 4.10.6, 

AG(kg) = -10 logjg g^ . (4.11.3) 

It is seen, from 4.9, 4.10, and 4.11, that the antenna parameters such 
as effective aperture area, gain, and gain degradation are defined using 
the total received power and not the dc power. Therefore, gain degrada- 
tion is not the same as the measured dc power degradation. The gain 

includes both fluctuating and non-fluctuating components of the incoming 
wave. Therefore, under conditions of heavy turbulence, it does not 
represent the usable power in the wave. 

It follows from the above discussion that gain degradation must 
represent the fluctuating power that is scattered out of the antenna 
beam. This will be discussed with examples in Chapter Seven. 

The two classes of receivers have now been delineated. A quick 

survey of the literature indicates that the distinction between the 

results for the two classes of receivers in the presence of turbulence 
has not been appreciated. Indeed, it appears that measurements of turbu- 
lence have been generally performed using asynchronous AM receivers, 
employing square law detectors. This clearly does not present a complete 
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picture of the turbulence as the phase information is lost. However, 
the more serious problem is that the statistics of the experimental 
results appear to be then calculated using the principles of the syn- 
chronous AM receiver! 

To summarize this chapter, two general classes of receivers, 
called synchronous and asynchronous receivers, respectively, have been 
defined. An important distinction has been drawn between the incoming 
signal, the output signal, and the measured signal or the incoming sig- 
nal as perceived by these receivers. Though the measured signal will 
be the same for both synchronous and asynchronous AM receivers in the 
absence of turbulence, fundamental differences occur when scintillation 
is present. Thus, care must be taken to ensure that the proper receiv- 
er is used for the intended application, and that the appropriate 
results are used when calculating signal characteristics in turbulence. 
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CHAPTER V 


THE CIRCULAR APERTURE 


5.1 General Relationships 

The generalized relationships for the problem of the aperture 
antenna in a turbulent atmosphere have been obtained in the preceding 
chapters. These results will be applied to the specific case of the 
circular aperture antenna in this chapter. Fundamental expressions will 
be‘ evaluated for any arbitrary angle of incidence of the incoming wave 
which does not deviate significantly from the main beam. The particular 
case of on-axis incidence will be then considered. The analysis will be 
carried out for a circular aperture with no feed taper, as well as for 
the Gaussian tapered aperture antenna. 

The circular aperture is of general interest. As shown in Chapter 
Two, the results of the formulation are independent of the mechanics of 
the feed system. Thus, the circular parabolic antenna can be studied, 
for example. Furthermore, with the proper choice of parameters, the 
results can be applied to the field in the focal plane of a focusing 
system. 

The geometry of the problem is shown in Figure 5.1. Any point on 
the plane aperture, diameter 'D', is specified by the cylindrical polar 
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z 



Figure 5.1. Geometry of the circular aperture antenna. 

coordinates (p,C). The radius 'a' is chosen to be the characteristic 
dimension. We define a normalized set of cylindrical coordinates, 

(R, 5) or where 

R = I . (5.1.1) 


Then, 


X = a R COS 5 
y = a R sine 

p 

The Jacobian of the transformation is a R. 
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The relative correlation length or correlation ratio is 


l 21 
P _ n i_n 

^ ~ T ~ D • 


(5.1.2) 


The distance between any two points ITj^ and R 2 is such that 


= a^(R^ + R^ - 2 RjR2COs(Cj - ^ 2 )), 


so that the correlation coefficient for a Gaussian model atmosphere is 


-(Rj^ + R 2 “ 2RjR2COs(5j - ^ 2 )) 


^ ® 


llj - r/ 


Also, let the generalized angle of arrival be 


(5.1.3) 


0 = -ka sine =-ir sine . 


(5.1.4) 


The vector (k ,k ), defined by (2.1.4), which indicates the angle 
A y 

of arrival, (0,4'), transforms to the corresponding generalized angle of 
arrival vector ( 0 , 41 ). 

The following are basic integrals which pertain to the circular 
aperture. From (4.2.3) 
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If this is evaluated by the power series method in (4.3.12), then 


1 2n 1 2 it 
0 0 0 0 


j0R,cos(C,-<^) -jeR^cosCCo-*^) 
e e ^ ^ R^R 2 dC^dR^d^ 2 dR . 


(5.1.10) 


For the Gaussian correlation coefficient, 


b 


m 

n 


-IRj-R2I^ 



9 


(5.1.11) 


where 



(5.1.12) 


If the magnitude of f(l^) is a function of | R| only and the direction of 
the polarization vector of the feed illumination, p^, is independent 
of position on the aperture, then all results will be independent of 4> 
and will be functions of the generalized angle of arrival, 0, only. 

From (5.1.10), (5.1.11), (4.3.12), and (4.3.13), it is seen that 

the only parameter that changes with m in Ip„(k ) is C . Therefore, we 

urn s m ’ 

define 
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' 'erne's) • 


(5.1.13) 


Ic(Cm,ks) win be called the correlation integral for the circular aper- 
ture. 


5.2 Uniform Illumination 


This section studies the uniformly illuminated circular aperture. 
The results are applicable to linear or circularly polarized focusing or 
non-focusing radio aperture antennas. A focusing antenna is one which, 
by using reflecting or refracting elements, directs the incoming rays to 
a focal point. The optical telescope, which may be considered to be a 
polarization insensitive circular focusing aperture, is another example 
of a uniform aperture. 

A uniformly illuminated aperture is represented by the aperture 
function 


f(R) = 1-pf 


(5.2.1) 


Then, from Appendix A5, the magnitude of the voltage gain pattern is 


9 ^( 0 ) = 


2Ji(0) 

0 ~ 


(5.2.2) 


Also, 


max 



(5.2.3) 
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(5.2.4) 


^2 ^ n 2 


^PO max ^ ^ ^ ^ ^ 


phys 


The plane wave power gain pattern function is 


9o(0) =1 


2Ji(0) 


'>2 




If 


6 =£® 

C 2 ’ 


the correlation integral is 


I (C 31 - VV 

p=0 q=0 


t=0 


Oc'c.®) ' Z 7^ 

\ 


r f - 

' r (^i)J 


where 


(5.2.5) 


(5.2.6) 


(5.2.7) 


(5.2.8) 
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V = p - 2q 
u = 2t + p + 1 


J ( 0 ) = the cylindrical Bessel function 

V 

S ( 0 ) = the Lommel function. 

5.3 On-Axis Incidence (uniform illumination) 

For on-axis incidence, 9=0. Hence, 

0 = 0 
9v«» " 1 
9o(0) = 1 


\ 


(5.3.1) 

(5.3.2) 



(5.3.3) 





\ 




m=0 


(5.3.4) 


It can be shown easily, by considering the individual terms, that 

Iq(C, 0) — > 1 as ® 

Ic(C,0)->C^ as C->0. 
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Figure 5.2 shows the axial incidence correlation integral, 
Iq(C, 0), as a function of the correlation ratio, C. The figure also 
shows an empirical approximation to Iq(C, 0) given by 


ic(c,o) ^(1 




1 

r 




(5.3.5) 


This expression was arrived at by inspection of the first few terms of 
(5.3.3) and by trial and error. It is a simpler expression to evaluate. 
However, the rigorous expression, (5.3.3), will be used in this study. 

It is seen that Iq acts like a high pass filter with a cut off 
point, at which its value falls to ^ being given by 


C 


crit 


■V 1 


(5.3.6) 


C . will be-called the critical relative correlation length. Thus, 

C I 1 u ^ 

when the correlation ratio is very high, gjj(C,a^,0)— > 1. From (4.3.8), 

When the correlation length, is comparable to the antenna 
size, (C 1), the gain degradation factor becomes progressively smaller. 
When the correlation length is much less than the antenna size, (C << 1). 

-o^ , _p2 

Then, g^(C,a^,0) e ^ and Pr-*Pro e 

Thus, the range of the received power is 



^'r - '^RO 


(5.3.7) 
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Figure 5.2. Correlation integral I^(C,0) for a uniformly illuminated 
circular receiving aperture antenna for on-axis incidence. 
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Figures 5.3 and 5.4 are used to illustrate these effects. Figure 5.3 
shows the behavior of on-axis gain, G(j(0) (4.10.6), at 30 GHz, as the 
antenna size is increased from 10 cm to 100 m. A correlation length, 

4^, of 10 m was assumed and the family of curves represents values of 
2 

o of 0, 1, 3, and 5, respectively. These numbers, while not unreason- 
able, are an example of a very turbulent atmosphere and, hence, have 
been chosen for illustrative purposes only. 

Figure 5.4 shows the gain degradation, aG( 0) (4.11.3), for the 
same set of conditions. It is, in fact, the behavior of the degradation 
integral, gd(C,o^,0). 


5.4 Gaussian Tapered Illumination 

The aperture distribution of the feed illumination is often 
tapered to obtain some desired property of the antenna. For example, 
where a low side lobe level is desired, the energy incident on the rim 
of the aperture, when considered as a transmitter, may be reduced to 
minimize edge diffraction effects. 


One distribution which is an example of this class of problems 
and which also lends itself to analytical studies is the Gaussian taper. 
The Gaussian taper circular aperture illumination function is represent- 
ed by 


R_ 

’t2. 

T(R) = e p. 


(5.4.1) 


Then, from Appendix A6, 


9v(0) ^ 



(5.4.2) 
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Figure 5.3. Gain of uniformly illuminated circular receiving aperture 
antennas in turbulence. 
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Figure 5.4. Gain degradation of uniformly illuminated circular receiv- 
ing aperture antennas in turbulence. 
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(5.4.3) 


max " ® ^ 


{' 


^PO max ~ 1 %hys ^ (1 “ ® ) 


(5.4.4) 


9o(0) = e ^ 


If 



(5.4.5) 


(5.4.6a) 


(5.4.6b) 


the correlation integral for Gaussian illumination is 


^(C,0) = 



Qc(Y,0)^ 

p! 


(5.4.7) 


Q0(y, 9) is defined in (5.2.8) with <5^ substituted for 
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5.5 On-Axis Incidence (Gaussian illumination) 


e = 0 


9v(0) = 1 

gQ(0) = 1 


(5.5.1) 

(5.5.2) 


Iy(C.O) 


, f./2\2pf2 , s|2 

= (5.5.3) 

.1-e^J P=0 


9d(^*V0) 


■°w V (‘'w)"' 


0 ) 


(5.5.4) 


m=0 


Figure 5.5 shows the correlation integral for Gaussian illumina- 
tion, Iy(C,0), as a function of the correlation ratio, C, for values of 
the taper, t, from 1.0 to 0.1. When t is large, the results tend to the 
uniform illumination case, as would be expected. Narrowing the aperture 
illumination function by reducing t slides the curve to the left, allow- 
ing the relative correlation, C, to be smaller before any degradation 
effects are noticeable. It will be seen from (5.4.2) that reducing t 
also widens the beamwidth. 

Reducing t, in essence, reduces the effective aperture of the 
antenna. The antenna appears to be smaller than its physical size alone 
would indicate. Therefore, for a given correlation ratio, C, the effec- 
tive correlation ratio, is larger than thus reducing degrada- 

tion effects. 
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Figure 5.5. Correlation integral I (C,0) for a Gaussian illuminated 

circular receiving aperture antenna for on-axis incidence. 
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Inspection of Figure 5.5 and Equation (5.4.6a) shows that if 
T is small, its effect would dominate over the effect of C, provided C 
is reasonably large. Also, Figures 5.5 and 5.2 are seen to be of the 
same form. Therefore, a simple relationship was sought between the 
uniform and Gaussian illumination cases. 

It was determined empirically, by inspection of Figure 5.5, that 
for T < 1, the critical value of C, i.e., that C at which the correla- 
tion integral was was given approximately by ■ 


Therefore, the equivalent uniform illumination correlation ratio is 


^eff ''' — 




T < 1 


(5.5.6) 


and, hence. 




'eff 




tD 


(5.5.7) 


This implies that the effective diameter 


Deff^V^D . (5.5.8) 

Thus, all results derived for the uniform illumination case may be 
used for the Gaussian illuminated aperture, provided that the effective 
diemeter, 0^^^, and correlation ratio, are used instead of D and C. 
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Finally, antenna designers usually prefer to work with the feed 
power taper, Tp, instead of the illumination taper constant x. Tp is 
defined as the ratio of the maximum value of the feed illumination 
power pattern, to the value in the direction of the rim of the aperture, 
and is expressed in dB. It is easily shown that. 


T 8.686 

-J- 


dB . (5.5.9) 


Summarizing this chapter, the general circular receiving aperture 
was evaluated for any arbitrary angle of arrival of a wave perturbed by 
a turbulent atmosphere, within the main beam of the antenna. Both 
uniform and Gaussian tapered feed illuminations were then completely 
evaluated. Next the special case of on-axis incidence was studied for 
both illuminations. The Gaussian tapered beam was related to the uni- 
form aperture by a simple expression. This permits the results obtained 
for the uniform aperture case to be applied to the Gaussian tapered 
circular aperture as well. 
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CHAPTER VI 


EXPERIMENTAL RESULTS 


In this chapter, the theoretical model developed in the previous 
chapters is matched to experimental work with satellite communications 
links. Experiments performed at The Ohio State University to measure 
the variance and the level of the received signal over different earth- 
space paths are used to assign numerical values to key parameters of the 
model. The model is then used to predict and is compared with the re- 
sults of similar measurements made elsewhere, to establish its validity 
under different conditions. 


6.1 General Remarks 


An important simplifying assumption used in this work is that 
the turbulence in the earth's atmosphere is homogeneous and isotropic 
over the path of the communications link. This is reasonable over short 
terrestrial paths.. Slant or earth-space paths encounter gradients in 
temperature, pressure, water vapor density and, hence, refractivity. It 
is necessary, therefore, to assume an equivalent homogeneous isotropic 
atmosphere for the purposes of the model, not only for earth-space paths, 
but also for comparison with terrestrial results. Its parameters would 
be, in some sense, an average of the parameters of the earth's atmosphere 
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over the path, for long term statistical purposes. The next two sections 
detail the experimental work performed at O.S.U. to determine these 
parameters. 

The primary experimental measurements used for this purpose were 
made at The Ohio State University ElectroScience Laboratory in Columbus, 
Ohio, when the Applications Technology Satellite, ATS-6, was moved from 
a geosynchronous position at 94° W. longitude to 35° E. longitude in 
1975, and back in 1976 [35,36]. The elevation angle of the satellite, 
viewed from Columbus, Ohio, (83° W. longitude) varied smoothly between 0° 
and 43°. The scintillation of the received signal amplitude as well as 
the average signal level were measured simultaneously at 30 GHz and 20 
GHz with a 4.5 m (15 feet) parabolic antenna during the descent (1975) 
and at 30 GHz (4.5 m antenna) and 2 GHz (9 m antenna) during the ascent 
(1976). The corresponding power taper, Tp, of the 4.5 m antenna was 
22 dB, while that of the 9.1 m antenna was 18 dB. 

The receivers employed square- law detectors and, hence, may be 
classed as asynchronous receivers. The voltage outputs were sampled at 
10 samples/second. The statistics of the equivalent antenna terminal 
voltage, v, as measured by the asynchronous AM receivers, were calculated 
for N samples after correcting for receiver calibration characteristics. 


6.2 Signal Amplitude Variance 

2 

The normalized experimental received signal amplitude variance, s , 
was calculated by dividing the amplitude variance by the dc power level 
and was expressed in dB below the dc power level. If v^- is the i*^" 
sample of v, from [ 35 ] 




(v^ - < v>)^ 
N <v>^ 


dB 


( 6 . 2 . 1 ) 


90 



where 


<v>=^^ . (6.2.2) 

i=l 

N = 2048 for the 1976 data. This corresponds to an averaging time of 
204.8 seconds for each set of samples. 

The data were taken under conditions of clear air at 33 elevation 
angles and the time average variances were calculated as the satellite 
ascended. The results were fitted to power law curves of the form 



A L' 


(6.2.3) 


where A and B are constants and L (km) is the equivalent path length 
through a homogeneous atmosphere. The atmospheric path lengths were 
calculated assuming a spherical earth of radius Rg, given by 4/3 the 
actual radius, or 8479 km, to correct for standard refraction, and an 
atmosphere of height, h. The equivalent path length is then given by 


L = 



+ 2hR^ + R^sin^e] - R^sine 
e e ■* e 


(6.2.4) 


where e is the elevation angle. 

A regressive fit was made to all the variance data to determine 
the values of A and B, while varying h to minimize the mean square error. 
The results at 30 and 2 GHz were [35], 
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2 in-6.4| 2.48 + .3 

30 - iu L 


(6.2.5) 


s2 ^ io-6-5l1.87± .2 


( 6 . 2 . 6 ) 


h = 6.0 + 1 km . (6.2.7) 

The fairly large error bounds do not indicate uncertainty in the data 
itself, but take into account the limited duration of the data periods 
and the correcting effect that observations over extended periods of time 
under different weather conditions would have on the results. 

The same analysis was performed on the 1975, ATS-6 descending, 
data [36], with N = 1024, i.e. , an averaging time of 102.4 seconds. The 
corresponding results for the 30 and 20 GHz data were. 


.2 _ in-6.2, 2.35 + .1 

^30 “ ^ 


h = 5.9 + 1 km . 


( 6 . 2 . 8 ) 

(6.2.9) 

( 6 . 2 . 10 ) 


Consequently, the height of the equivalent homogeneous atmosphere in the 

present model was also taken as 6 km. The signal amplitude variance as 

calculated by the present model (4.7.2) was then fitted in the minimum 

mean-square-error sense to the power law curves (6.2.5) - (6.2.9). 

The correlation length of the atmosphere, £ , and the variance of the 

2 ^ 

refractive index fluctuations, Op, were varied simultaneously to minimize 
the deviation from the power law curves. Subsequently, the height of the 
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equivalent atmosphere, h, for the model was also varied in an effort to 
further minimize the deviation. However, it was found that 6 km was the 
best overall compromise. The final results were 


0 ^ = 4.0 X 10'^^ 


= 46 m 
h = 6.0 km . 


( 6 . 2 . 11 ) 

( 6 . 2 . 12 ) 

(6.2.13) 


These results will be used in the subsequent discussion unless stated 
otherwise. 

The question now arises as to the reasonableness of these values. 

A search through the literature provided very few direct statements of 

p 

these values. Shifrin [37] indicates that o“ could take values from 

2.5 X 10"^^ to 2.5 X 10"^^. Ishimaru[28] indicates that is of the 

- 1 ? ” 
order of 10 . Shifrin also gives the value of 2.^ as ranging from 5 m 

to 100 m. 

Some work has been done, however, to measure the atmospheric struc- 

p 

ture constant, C^, of the Kolmogorov spectrum. Chadwick and Moran [ 38] 
and Gossard [39,40] are examples of recent publications. Height distribu- 
tions are also given in [39] . An alternative method used, therefore, 

p 

was to calculate the value of the equivalent structure constant, for 
the Gaussian spectrum, as defined in (3.7.3) and compare this with the 
published measurements instead. Then, from (6.2.11), (6.2.12) and 
(3.7.3), 

_2 

= 5.3 X 10'^^ m ^ . (6.2.13) 
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It is commonly accepted that near the surface varies over about 

n 

two orders of magnitude. Average values tend to vary between and 

depending on the height at which the measurement is made. The 
above value of the equivalent structure constant appears to be reason- 
able, especially since it represents an average over a wide range of 
path lengths and heights. 

Figure 6.1 shows the predicted signal amplitude variance at 2, 20 
and 30 GHz using the appropriate antennas, together with the power law 
approximations to the 1975 and 1976 data. The agreement is reasonably 
good over a range of frequencies of 15 to 1 and is a further confirma- 
tion of the values of the chosen parameters. 

It must be remembered, however, that the power law curves are 
approximations only and express only a path length dependence. They 
cannot, therefore, fully represent the variation of the signal variance 
with elevation angle. Figures 6.2 and 6.3 show the actual variances 
measured during the 1976 experiment, at each elevation angle, at 30 GHz 
and 2 GHz, respectively. The maximum and minimum observed values are 
plotted, together with the long term average predicted by the model. 

The results show a good fit. Indeed, the agreement here appears to be 
better than that shown in Figure 6.1. 

As a further check, the model was used to calculate the variance 
for the IDCSP X-band beacon measurements reported by R.K. Crane [12]. A 
60-foot parabolic antenna was used at 7.3 GHz over elevation angles 
from 0.5 to 10 degrees. The results were recast to be consistent with 
the definition of variance in this report (D.M. Theobold [36]). The 
prediction is shown in Figure 6.4. The agreement is again quite good. 

Other measurements were made at The Ohio State University Electro- 
Science Laboratory using the Communications Technology Satellite, CTS, 
at 11.7 GHz and the COMSTAR Satellite at 28.56 GHz [41]. Measured 
and predicted amplitude variances are summarized in Table 1. The CTS 
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Figure 6.1. ATS-6 measured amplitude variance at 2, 20, and 30 GHz com- 
pared to theoretical model. 
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Figure 6.2. ATS-6/1976 limits of measured amplitude variance at 2 GHz 

compared to theoretical model. 
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Figure 6.3. ATS-6/1976 limits of measured amplitude variance at 30 GHz 

compared to theoretical model. 
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Figure 6.4. IDCSP X-bandmean and limits of measured amplitude variance 
compared to theoretical model. 
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TABLE 1 

SOME COMPARISONS OF THE THEORY AND EXPERIMENTAL RESULTS 


Satellite 

Frequency 

(GHz) 

Receiver 
Antenna 
diam (m) 

Elevation 

Degrees 

Variance 

(dB) 

Measured Model 

CTS 

11.7 

0.6 

32.7 

- 36 

- 42 

COMSTAR D 3 

28.56 

0.6 

43 

- 42 

- 45 


measurements were heavily contaminated by system noise due to the low 
link margin. The other measurement with the 0.6 m antenna is also 
influenced by system noise to some degree. Under the circumstances, the 
agreement is reasonable. 

Thus, the model appears to predict the signal amplitude variance, 
s^, satisfactorily as a function of frequency, path length, and aperture 
size. 


6.3 Signal Level Degradation 

The experimental received dc signal power, Pj^, as measured by the 
receiver was defined as 

= <v>^ (6.3.1) 

where < v> has been defined in (6.2.2). Since the receivers used in these 
experiments incorporated square-law detectors, the corresponding 
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theoretical result for the received power level is defined by (4.6.3). 


The average received dc power was computed as a function of eleva- 
tion angle from (4.6.3) for the receiving systems used in the ATS-6 

I 

return measurements [35]. Gas loss was added to account for the attenua- 
tion caused by the atmosphere. The 6 km equivalent atmosphere was 
assumed for this purpose. The results were normalized to the power 
that would be received at zenith and are shown in Figure 6.5. 

In addition to the 2 GHz and 30 GHz data, median signal levels as 
a function of elevation angle, measured by McCormick and Maynard at the 
Communications Research Center, in Ottawa, using the US TACSATCOM-1 
7.3 GHz beacon [42] are also available. A 9.1 m antenna was used. 
Measurements were made over elevation angles from 6° to 0.5° as the 
satellite moved eastward. The data has been recast to agree with the 
definitions used in this report by Theobold [36]. These measurements 
and predicted results are also included in Figure 6.5. 

The theoretical results were calculated at 2, 7.3 and 30 GHz for 

9.1, 9.1, and 4.5 m antennas, respectively. These correspond to 60, 

220, and 450 wavelength antennas, with relative correlation ratios of 

10.1, 10.1, and 20.4, respectively. 

The agreement is seen to be quite good. The model thus appears 
to predict long term average power level degradation as well. 

Figure 6.6. shows the average power degradation, as calculated 
from (4.8.2) only, excluding the gas loss. When compared with Figure 
6.5, gas loss is seen to be the dominant cause of reduction in power 
levels at low elevation angles. However, the turbulence contribution is 
also significant. It will also be shown later that the effect can be 
considerably larger in millimeter wave systems. 
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Figure 6.5. ATS-6 and TACSATCOM-1 measured received signal levels com- 
pared to theoretical model. 
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Figure 6.6. ATS-6 and TACSATCOM-1 degradation component of received sig- 
nal level. 
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In summary, numerical values were derived for the constants of the 
model by fitting it to available experimental data. The data used were 
from different locations. The relevant parameters are given by (6.2.11)- 
(6.2.13). The model predicts long term average signal variance and 
received power over a wide range of frequencies and antenna sizes. 
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CHAPTER VII 


DISCUSSION OF THE MODEL 


7.1 Introduction 


The key parameters of the long term average behavior of the model 
have now been estimated. Using these estimates, cases of interest are 
studied to gain insight into the behavior of the model as well as to 
predict the properties of the communications system. 

Four frequencies have been chosen for this purpose. These are 3, 
12, 30 and 90 GHz. They span the range of current interest and, in 
the case of 90 GHz, represent a frequency band in which increasing in- 
terest is being shown. While the validity of the model at 90 GHz is, at 
present, unknown it should provide at least an initial estimate for the 
link designer. As data become available, they could be easily compared 
with the given curves for verification. 

Four representative antenna sizes have also been chosen for de- 
tailed study. These are 1 m (3.3 ft.), 4.5 m (15 ft.), 9.1 m (30 ft.), 
and 30 m (100 ft.), respectively. The first three are popular sizes. 

The 30 m antenna has been used at the lower microwave frequencies, es- 
pecially in class A earth stations. It is unlikely that this size will 
be used extensively at 90 GHz, as the required tolerances become pro- 
hibitive. However, it represents an extreme case and will be a severe 
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test of the soundness of the model. It will also allow some properties 
of the propagation link, which might otherwise be hidden, to be brought 
out. All results presented at 90 GHz for the 30 m antenna should be 
viewed in this light. Also, whenever possible, the parameter under 
study is plotted at 30 GHz with the antenna aperture size varying con- 
tinuously from 10 cm to 100 m. Obviously, the comments made above apply 
to the larger sizes in these plots as well. 

Average conditions of turbulence will be always assumed, unless 
specifically stated otherwise. However, some properties of the link 
such as antenna gain degradation, which may not be significant under 
quiet conditions, may have a major impact on the performance of the 
system under conditions of strong turbulence. Hence, some cases are 
evaluated with this in mind, in order to estimate the worst case behav- 
ior of the link. 

Some general properties of the formulated model atmosphere will 
be first discussed. Next, amplitude variance will be studied, both for 
the synchronous and the asynchronous AM receiver. The study will be 
repeated for the degradation of the average received signal power. 
Finally, the gain degradation of the antenna will be examined under 
conditions of strong turbulence. 


7.2 The Atmosphere 


The long term average parameters of the atmosphere, modeled as an 
equivalent homogeneous isotropic Gaussian spectrum, were estimated to be, 
(6.2.11) - (6.2.13), 

al = 0.4 X 10“12 (7.2.1) 

Af, = 46 m (7.2.2) 
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h = 6 m 


(7.2.3) 


The corresponding correlation coefficient is shown in Figure 7.1. 
The figure also shows the correlation function of the equivalent Von 
Karman spectrum for homogeneous turbulence, having the same correlation 
length. Figure 7.2 shows the corresponding spectra for the two models. 

If the Von Karman spectrum is accepted as being closer to the natural 
spectrum, it will be seen that errors may occur in the Gaussian model in 
the geometrical optics region and in the long path region of a propagat- 
ing wave. 

Figure 7.3 shows the path length through this 6 km equivalent 
atmosphere over a 4/3 radius earth, for elevation angles from 0° to 45°. 
The maximum length (at 0°) is about 320 km. Thus, the ratio of the path 
length along the horizon to that at zenith is, approximately, 53. There- 
fore, a power quantity which is directly proportional to the path length 
will vary by about 17 dB from zenith to the horizon due to the variation 
in path length alone. 

Figure 7.4 shows the wave variances (3.5.13b) of waves at the four 

frequencies of interest, namely 3, 12, 30, and 90 GHz, as they propagate 

through this atmosphere. It should be remembered that the wave variance 

is an artifical construct of no physical significance, being defined 

simply as the sum of the log amplitude and phase variances. The path 

length characteristic of Figure 7.3, is clearly mirrored in these curves, 

2 

as expected. The k dependence of values is also seen. 

Figure 7.5, showing the phase variance of the propagating waves, 
is very similar to Figure 7.4. This is because the long path region has 
not been reached. Consequently, the phase variance is much larger than 
the log-amplitude variance and dominates in the wave variance. Some 
differences are becoming noticeable at 3 GHz at low elevation angles. 


106 



CORRELRTION COEFFICIENT 
0.01 0.10 1.00 



SEPfiRfiTION (M) 


Figure 7.1. Correlation coefficient of the model atmosphere and the 
equivalent Von Karman model. 

Op = 0.4x10“^^; Ap = 46 m; Lq = 64 m; Hg = 10 mm. 
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Figure 7.2. Spectrum of the model atmosphere and the equivalent Von 
Karman model. 

= 0.4x10”^^; = 46 m; Lq = 64 m; Hg = 10 mm. 


108 




jh the model 


6 km atmosphere. 


9 




WAVE VflRIRNCE (DB) 

-60 -50 -40 -30 -20 -]0 


PftTH LENGTH (KM) 

320 100 50 3C 20 ]5 12 10 9 



ELEVfiTION (DEG) 


Figure 7.4. Wave variance dependence on frequency. 
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Figure 7.5. Wave phase variance dependence on 
frequency. 
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The log-amplitude variance is shown in Figure 7.6. The character- 
istics are markedly different from those of the previous two variances. 
Numerical values are several orders of magnitude below the phase vari- 
ances. Thus, it may be expected that the phase ripple across the wave 
front should be a dominating feature of the wave. Consequently, there 
would be significant phase perturbation of the antenna terminal voltage 
as well as some amplitude scintillation caused by the incoming rays not 
adding in phase. Both of these effects would affect the scintillation 
measured by a synchronous, or phase sensitive, receiver. 

A receiver whose demodulating circuits are sensitive to amplitude 
scintillations alone should, hence, perceive a much smaller scintilla- 
tion. This explains the curves of measured scintillation in Chapter Six 


7 , 3 Signal Amplit u de Variance 

Figure 7.7 shows the variance of the amplitude scintillations 
measured by an asynchronous receiver at the four frequencies of interest 
The diameter of the circular aperture has been fixed at 4.5 m for these 
curves. It is interesting to note that the variance is beginning to 
saturate at low elevation angles. This will be discussed in greater 
detail shortly. 

Figure 7.8 shows the corresponding curves for a, system with a 
synchronous receiver. The variance seen by this receiver is several 
orders of magnitude higher than the preceding case for the reasons 
explained before, i.e., the phase perturbations across the aperture 
introduce amplitude scintillation and phase fluctuation of the antenna 
terminal voltage. Therefore, the synchronous receiver, which is sensi- 
tive to phase fluctuations as well, would see a much higher fluctuation 
of the terminal voltage. The characteristics of the phase fluctuations 
clearly dominate the results. 
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Figure 7.6. Wave amplitude variance dependence on frequency. 
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Figure 7.8. Frequency dependence of signal amplitude variance, measured 
by a synchronous receiver. D=4.5 m, F=3, 12, 30 and 90 GHz. 
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Curves for the variance of the voltage measured by synchronous 
receivers also may be seen to show the variance of the field in the 
focal plane of a focusing antenna, as noted in Section (2.6). It should 
be useful to bear this in mind during this discussion. 

The results shown above are for long term average conditions only. 
An attempt was made to estimate results for extreme conditions. It was 
felt that the extreme cases in the atmosphere occur first under 
very quiet conditions, i.e., when the refractive index variance is at 
a minimum. It would then appear, intuitively, that the correlation 
length would be very high. The other extreme would be when is very 
high. It is reasonable to expect that the correlation length should 
then be low. This would also pose the worst case for the antenna. To 
obtain these conditions, the refractive index variance was changed by a 
factor of 10 about its median value and the correlation length was also 
varied by a factor of 10. The corresponding values of were also 
computed for comparison purposes. 


Case 1: Quiet atmosphere 


0 ^ = .04 X 10"^^ 


= 100 m 

<^ng ■ -3 * 



Case 2: Disturbed atmosphere 


= 4.0 X 10"^^ 
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= 10 m 

= 1.4 X 10"^^ 

ng 
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The equivalent structure constant is for comparison purposes only. 
The results are shown for the two receiver systems in Figures 7.9 and 
7.10 at 30 GHz for an aperature diameter of 4.5 m. These may be reason- 
ably expected to be the extreme conditions to which the links could be 
expected to be subjected under clear air conditions. 

The validity of the present model under conditions of strong tur- 
bulence is not known. Much depends on the validity under these condi- 
tions of the amplitude and phase filter functions used to calculate the 
log amplitude and phase variances. Furthermore, the expressions for the 
average power and signal variance measured by an asynchronous receiver 
may not contain a sufficient number of terms of the Taylor series expan- 
sion used to remain accurate. An injudicious choice of the refractive 
index variance and correlation length can cause these expressions to 
blow up. Unfortunately, improving their accuracy necessitates the 
evaluation of ever higher moments of the asynchronous received voltage. 

A reference to Appendix A4 quickly discourages this approach. 

Therefore, all results shown for the strong turbulence case must 
be used with care until their validity is verified by theoretical or 
experimental methods. 

The dependence of amplitude variance on aperature size is examined 
in Figure 7.11 for the asynchronous receiver case. Only 1 m and 9.1 m 
diameter apertures, which encompass the range of popular sizes, are 
considered at each frequency. The results are shown separately in Fig- 
ures 7.12 to 7.15. There is a clear increase in the measured scintilla- 
tion with the larger aperture. 
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Figure 7.9. 


Change of variance with extremes of turbulence, F=30 GHz; 
D=4.5 m; asynchronous receiver. 

2 -12 

Average turbulence: = 0.4x10 ; = 46 m. 

2 -12 

Minimum turbulence: = .04x10 ; = 100 m. 

Maximum turbulence: = 4.0x10“^^; = 10 m. 
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Figure 7.10. Change of variance with extremes of turbulence, F=30 GHz; 
0=4.5 m; synchronous receiver. 
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Figure 7.11. Dependence of variance on antenna aperture size, at 3, 12, 
30 and 90 GHz; asynchronous receiver. 
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Figure 7.12. Dependence of variance on antenna aperture size, at 3 GHz; 
asynchronous receiver. 
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Figure 7.14. Dependence of variance on antenna aperture size, at 30 GHz 
asynchronous receiver. 
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Figure 7.15. Dependence of variance on antenna aperture size, at 90 GHz; 
asynchronous receiver. 
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Figure 7.16 is a different presentation of the variance of the 
received signal measured by an asynchronous AM receiver. The measured 
amplitude variance is displayed in the figure as a function of aperture 
size, at selected elevation angles. The path length (km) through the 
equivalent homogeneous atmosphere is given in parentheses. The atmos- 
phere is assumed to be under conditions of average turbulence. 

The normalized amplitude variance is initially constant at any 
given elevation angle. However, it increases rapidly with the size of 
the aperture beyond some critical value in the range of 1 m to 10 m. 
Another important characteristic is the evidence of aperture averaging 
at the larger aperture sizes over the longer propagation paths. There 
is a clear reduction in the increase of the amplitude variance with 
aperture size, compared with that at higher elevation angles. 

It is essential to note that this should not be classed simply 
as a D/X effect. A fundamental feature of the present model is that the 
diameter of the antenna is related not only to the wave length, X, but 
also to the size of the correlation length of the atmosphere. Consequent- 
ly, two systems, having the same electrical aperture size, but at differ- 
ent frequencies, may have marked differences in behavior because of 
the different physical aperture sizes and, hence, different relative 
correlation ratios to the atmospheric turbulence scale size. In other 
words, since the incoming wave front is no longer plane, the number of 
ripples across the aperture would be different at different frequencies 
for antennas of the same electrical aperture size, thus affecting the 
results by different amounts. 

The asynchronous receiver eliminates the phase information in the 
measured signal. Therefore, it is difficult to relate the observed 
results directly to the properties of the received wave. Equation (4.7.2) 
exemplifies this problem. The synchronous receiver, on the other hand, 
retains more information from the received signal. Hence, the amplitude 
variance measured by a synchronous AM receiver is studied next. 
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Figure 7.16. Dependence of variance on antenna aperture size, at select- 
ed elevation angles; F=30 GFIz; asynchronous receiver. The 
equivalent path length is shown in parentheses. 
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Figure 7.17 shows the amplitude variance as a function of elevation 
angle for aperture sizes of 1 m and 9.1 m only, at the four selected 
frequencies for synchronous receiver systems. Average conditions of 
turbulence are assumed as before. There is no noticeable effect on the 
variance with change in aperture sizes in this range. Therefore, the 
normalized amplitude variance (4.4.3) was next plotted in Figure 7.18 for 
aperture sizes from 10 cm to 100 m at 30 GHz, as in Figure 7.16. 

While there is the expected change of the variance with elevation 
angle, the aperture size dependence is weaker. Furthermore, instead of 
increasing with aperture size, as in Figure 7.16, the variance remains 
constant for small apertures and then actually decreases with larger aper- 
ture sizes. 

Unfortunately, the normalization of the variance to the measured 
dc power, while useful as an experimental procedure, introduces an addi- 
tional variable, thereby making comparison of Figures 7.16 and 7.18 and 
the interpretation of the results difficult. The true fluctuating power 
in the received signal, ^fl (4.3.19), is a more useful aid to understand- 
ing the mechanisms involved. 

Figure 7.19 shows the measured fluctuating power, P^-|, normalized 
to the plane wave power, received in the absence of turbulence. The 
latter is, by definition, independent of the turbulence and the elevation 
angle. P^-| cannot, obviously, exceed Pj^g. The curves are drawn under 
conditions of average turbulence in the same format as Figures 7.16 and 
7.18. 


At high elevation angles the fluctuating power increases with path 
length as expected. However, an interesting behavior is revealed by the 
0° and 1° curves, which actually cross over and indicate smaller values 
of the received fluctuating power at long path lengths for a given large 
aperture size, than at shorter path lengths. To bring out this effect, the 
same cases are plotted for maximum turbulence conditions in Figure 7.20, 
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Figure 7.17. Dependence of variance on antenna aperture size, at 3, 12, 
30 and 90 GHz; synchronous receiver. 
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Figure 7.18. Dependence of variance on antenna aperture size, at select- 
ed elevation angles; F=30 GHz; synchronous receiver, the 
equivalent path length is shown in parentheses. 
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Figure 7.19. Dependence of the fluctuating component of the received 
power on antenna aperture size at selected elevation 
angles in average turbulence; F=30 GHz; synchronous re- 
ceiver. 

= 0.4x 10'12. = 45 m . 

The equivalent path length is shown in parentheses. 
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Figure 7.20. Dependence of the fluctuating component of the received 
power, on antenna aperture size at selected elevation 
angles in maximum turbulence; F=30 GHz; synchronous re- 
ceiver. 

= 4.0x10’^^; = 10 m. 

The equivalent path length is shown in parentheses. 
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where the phenomenon is magnified and visible at shorter path lengths as 
well. It is also apparent that the crossover occurs at antenna sizes 
close to the size of . the correlation length of the refractive index fluc- 
tuations. It is suggested that the path length at which the effect 
becomes noticeable is such that the scale size of the turbulence becomes 
comparable to the size of the Fresnel zone of the antenna. 

Figures 7.16 through 7.20 illustrate the phenomenon known as aper- 
ture averaging. Over short path lengths, the phase fluctuations are well 
correlated over the Fresnel zone of the antenna. Therefore, increasing 
the path length, which increases the magnitude of the amplitude and phase 
fluctuations results in an increased signal variance. However, when the 
path length is sufficiently large and the antenna is also large, the 
phase correlation length becomes comparable to or less than the Fresnel 
size of the antenna. Now, different rays reaching the antenna are uncorre- 
lated, causing partial cancellation of the fluctuating component of the 
signal. This also reduces the variance of the antenna terminal voltage. 

An alternative viewpoint is to look at the antenna as a spatial 
filter. Referring to Chapter 4.3 - 4.4, it is seen that if the antenna 
is small compared to the scale size of the turbulence, the amplitude 
variance measured by a synchronous receiver is, in fact, the true fluctu- 
ating power in the wave. This, clearly, will be a function of the atmos- 
pheric turbulence only, as long as there is no significant scattering out 
of the beam of the antenna. 
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1 . 4 Signal Level Degradation 


The degradation in received power as perceived by the receiver 
was examined under the same conditions and for the same dependences as 
the amplitude variance. The following is a summary of the findings of 
the resulting degradation. 

Figures 7.21 and 7.22 show the frequency dependence of the power 
degradation for the asynchronous and synchronous receiving systems, re- 
spectively. The former does not show significant effects except at very 
low elevation angles. The latter, as shown before, measures the true 
non-fluctuating component of the incoming wave. Thus, degradation ef- 
fects are significant at all elevation angles. 

The aperture size dependence of the average power degradation 
experienced by an asynchronous receiver at 30 GHz is shown in Figures 
7.23 and 7.24. It is displayed as a function of elevation angle in 
Figure 7.23 for the maximum turbulence case. It is also plotted with 
aperture size as the primary variable in Figure 7.24 for the average 
turbulence case. Even if the 0° values are excluded, it is significant 
that there may be more than a tenfold decrease in the received average 
power at low elevation angles in strong turbulence, over that received 
under average turbulence conditions. Hence, system margins may easily be 
degraded in heavy turbulence. The likelihood of such an occurrence is, 
of course, not estimated in. this report. 

Average power degradation for a synchronous receiver is shown in 
Figure 7.25. It is included for completeness and for purposes of compar- 
ison only. There is no aperture size dependence of power degradation for 
synchronous receiver systems. This is seen in Equation (4.3.17). The dc 
power measured by the synchronous receiver is the true coherent power 
in the wave. This is, clearly, independent of the receiving system and 
is a function of the atmospheric link only. 
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Figure 7.21. Frequency dependence of average power degradation, D=4.5 m; 
asynchronous receiver. 
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Figure 7.22. Frequency dependence of average power degradation, D=4.5m 
synchronous receiver. 
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Figure 7.23. 


Aperature size dependence of average power degradation, in 
maximum turbulence; F=30 GHz; asynchronous receiver. 

al = 4.0x10“^^; = 10 m. 
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Figure 7.24. Dependence of the average received power on antenna aper- 
ture size at selected elevation angles in average 
turbulence, F=30 GHz; asynchronous receiver. The equival- 
ent path length is shown in parentheses. 
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Figure 7.25. Dependence of the average received power on antenna aper- 
ture size at selected elevation angles in average 
turbulence, F=30 GHz; synchronous receiver. The equivalent 
path length is shown in parentheses. 
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7.5 Gain Degradation 




The gain degradation, defined in Chapter 4.11, is examined briefly 
in this section under conditions of strong turbulence. Gain degradation 
has also been described as an 'aperture to medium coupling loss'. The 
gain is a measure of the total received power, i.e., the average as well 

as the fluctuating power is included in its definition. Gain degradation, 
hence, implies that some of the fluctuating power is not being added back 
into the output of the antenna. This can occur when the magnitude of the 

turbulence is such that some of the fluctuating power is scattered out of 

the beam of the antenna. Hence, a narrow beam antenna will suffer a 
greater gain degradation for a given degree of turbulence. 

In contrast, it should be recalled that signal level degradation 

was defined as the reduction in only the dc component of the measured 

signal. It is affected by the type of receiver used as well. Gain de- 
gradation, however, is independent of the receiver system and depends on 
the size of the antenna and the characteristics of the turbulence only. 

The frequency dependence of gain degradation is shown in Figure 
7.26 with elevation angle, under conditions of strong turbulence. 

The aperture diameter is 4.5 meters. Figure 7.27 shows the aperture size 
dependence at 30 GHz, also under conditions of strong turbulence. 

The aperture size dependence is displayed in a different format in 
Figure 7.28. The antenna gain is shown at 30 GHz as a function of eleva- 
tion angle (or equivalent path length). A region of gain saturation is 
seen. Increasing the aperture size in this region does not result in an 
increase in gain. However, a further increase in size narrows the beam 
so much that almost none of the fluctuating power is captured by the 
antenna. Then the gain resumes its increase with aperture size, now 
receiving the coherent power in the wave only. 
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Figure 7.26. Frequency dependence of antenna gain degradation, in maxi- 
mum turbulence. 

D = 4.5 m; = 4.0x10”^^; = 10 m. 
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Figure 7.27. Aperture size dependence of antenna gain degradation in 
maximum turbulence. 


F=30 GHz; = 4.0x10'^^; = 10 m. 


141 




Figure 7.28. Dependence of antenna gain on aperture size, at selected 
elevation angles in maximum turbulence. 

F=30 GHz; = 4.0x10"^^; = 10 m. 

The equivalent path length is shown in parentheses. 
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Figure 7.29 shows the same case as Figure 7.28, except that only the 
degradation is shown, instead of the total gain. Gain degradation in 
excess of 10 dB is possible, depending on the antenna size. 

It must be reiterated that the results shown in this chapter should 
be used with caution. The Gaussian approximation to the Von Karman spec- 
trum carries inherent inaccuracies. The expressions used may not be 
completely valid under conditions of strong turbulence. The results 
should be viewed primarily as illustrations of the processes taking place. 
However, the model permits at least an order of magnitude estimation of 
the effects of turbulence on the communications link. Its performance 
would be considerably better for predicting long term average behavior. 

In general, the use of numerical values should be combined with normal 
engineering practice, appropriate safety factors, and a healthy dose of 
skepticism! 
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Figure 7.29. Dependence of antenna gain degradation on aperture size, 
at selected elevation angles, in maximum turbulence. 

F=30 GHz; = 4.0x10"^^; = 10 m. 

The equivalent path length is shown in parentheses. 
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CHAPTER VIII 


SUMMARY 


The effects of atmospheric turbulence on microwave and millimeter 
wave communications systems, with special emphasis on satellite communi- 
cations systems, has been examined. 

A generalized theory of the receiving aperture antenna was first 
formulated. The results were extended to the case of an incident wave 
propagating through a turbulent medium. 

The earth's atmosphere was considered next. The Von Karman model 
was discussed and approximated by the Gaussian spectrum. The behavior 
of the log amplitude and phase of a wave propagating through such an 
atmosphere were examined. Some new results were formulated for the 
correlation functions and variances of these statistical parameters. 

The role of the receiver was studied next. It was shown that 
the effect of the receiver cannot be ignored in the presence of turbu- 
lence. The statistics of the received signal may be perceived very 
differently by the two general classes of receivers described as synchron 
ous and asynchronous receivers, respectively. The total power, variance. 
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signal level and the degradation ofthe signal level measured by each of 
these receivers were then formulated. It was noted that the results for 
the synchronous receiver may also be used to study the field at the focal 
plane of a focusing antenna. 

The particular case of the circular aperture was then considered. 

A complete result for the power received and the antenna gain was derived 
for both the uniform feed function and Gaussian tapered feed functions, 
provided the angle of arrival of the wave did not deviate significantly 
from the main beam. The particular case of on-axis incidence was studied 
in detail and some general characteristics of gain and gain degradation 
in turbulence were discussed. The Gaussian tapered aperture and the 
uniform aperture were shown to be simply related if the feed taper in the 
former is reasonably sharp. 

Experimental results were next used to determine the numerical 

values of key parameters in the atmospheric model. It was shown that the 

atmosphere could be approximated by a uniform homogeneous isotropically 

turbulent atmosphere over the propagation path. Its equivalent height 

was found to be 6 km over an earth of effective radius equal to 4/3 the 

physical radius, or 8479 km. The variance of refractive index fluctua- 
-12 

tions of 0.4 X 10 , together with a correlation length of 46 meters, 

was found to be suitable for prediction of long term average behavior. 

The model adequately predicted link behavior over a frequency ratio of 15 
to 1 and in different locations. 

The link model was used to study signal variance, degradation 
of average power, and antenna gain degradation for frequencies from 
3 GHz to 90 GHz and aperture diameters from 1 meter to 30 meters. The 
average power degradation indicates the reduction in the time average 
received signal level as measured by the receiving system. The gain 
degradation, on the other hand, is a measure of the reduction in the 
total power captured by the antenna. Their behavior was compared at the 
extremes of intensity of turbulence as well as for long term average 
values. 
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In general, the received signal variance increases with frequency. 
Its frequency dependence varies with the type of receiver used. The 
received power degradation also increases with frequency and changes with 
the type of receiver used. There is, however, no aperture size depend- 
ence of received power degradation when a synchronous receiver is used. 

Regions of gain saturation are seen with large antennas. Examples 
of aperture averaging were also found when the antenna size and the 
Fresnel size of the propagation path were comparable to the scale size of 
the refractive index fluctuations. 

It was shown that the synchronous receiver indicates a greater 
level of signal variance than the asynchronous receiver. Thus, the 
latter is preferable for communications links. However, the former 
should be preferably used to study atmospheric turbulence as its output 
contains phase as well as amplitude information. Care should also be 
taken to use the expressions appropriate to the particular receiver 
for evaluating signal level and variance. 

Finally, the numerical results should be used with care and a clear 
understanding of the approximations involved. Normal engineering prac- 
tice should be followed when applying the contents to system design. 
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APPENDIX A1 
MODEL SUMMARY 


Al.l Field Quantities 

The antenna aperture, S, lies in the x-y plane. P” = (x,y) is any 
point on this surface. A plane wave incident on the turbulent atmosphere 
and travelling a distance L through it, reaches the aperture at a direc- 
tion ( 6, ({>) with respect to the positive z axis. 

Let 

kg = (k^.ky) , (Al.1.1) 

where 


k^ = k sine cos<l> (Al.l. 2) 
ky = k sine sin<^ . (Al.l. 3) 
The incident field produces a random aperture field distribution 


?r(P) = S(F) e 



•P 


U(P) = Aq 


X(P)+j'l>(P) 


e 


(Al.l. 4) 


(Al.l. 5) 
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(Al.1.6) 


U(P) = p U(P) 

X = log amplitude 
^ = phase 


Let 


U„ = 


Aq e 




X„ = < X > 


(Al.1.7) 

(Al.1.8) 


The aperture illumination is described by the feed function f(P), 


f(P) = Pf f{P) 


(Al.1.9) 


Both the polarization of the incident wave, p, and the illumination, p^, 
are assumed to be independent of position at all times. 


A1.2 The Atmosphere 

The turbulence is homogeneous and isotropic over the path. The 
correlation function of the refractive index fluctuations is 





(Al.2.1) 
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/Vp 

C = ^ = The correlation ratio, is the correlation length of the 
^ turbulence, il^, normalized to a characteristic antenna 
dimension, L^. 


Also, 


4 = ^ 


tan"^W 


1 - 


W, 


2 _ YT 2. .2, 

o^ !5- o„i,„k L 


■^n~n 


tan 

1 + 9 

^ W 

9 J 


<^w = 


2,2 
O + 0 , 

X ^ 




o„2(l - b„) - 0^(1 - b^) + <,2(1 - b^) 


(Al.2.2) 

(Al.2.3) 

(Al.2.4) 

(Al.2.5) 

(Al.2.6) 


W 

g 



(Al.2.7) 


where 

2 

= variance of atmospheric refracted index fluctuations. 

p 

= log amplitude variance. 

2 

= phase variance. 
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= wave variance, 
w 


Wg = Gaussian wave parameter. 


b = log amplitude correlation coefficient. 
X 


= phase correlation coefficient. 


Average Received Voltage 


A1.3 Synchronous Receiver 


< V3y, > = < V,(k,) > 


1 2 

“^CFw - 

Vdo ® P*Pf . 


(Al.3.1) 


where 


9v(ks) 


ivC<s) 


I. 


V max 

. ' i 

is the voltage pattern function. 


(Al.3.2) 


“ j k * P 

Iv(ks^ = J n?) e ^ dS 

S 


(Al.3.3) 


''dO " max P*Pf 


(Al.3.4) 
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Asynchronous Receiver 


^ 'X/ 

< V > = < Iv He ) I > 

asyn I 


4a" 


IVdollgvC^s)! + 9d]4[^ ■ 29d + e ^]\ Ip.p^l 


(Al.3.5) 


is defined in (Al.4.3). It is the gain degradation factor. 


A1.4 Total Received Power 


~ ^RO ^0^*^s^^d |P*Pfl * (Al.4.1) 

where 

gpdj) = |gvffs)|^ (Al.4.2) 


is the power pattern function. 


9d 


-% w 


Ipo(ks) 


(Al.4.3) 


IpCk^) = 



/MW 


dSjdS2 


S S 


(Al.4.4) 
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Ipo(l<s) - 


(Al.4.5) 



_ * -jk ’(Pi-Po) 
f(Pl)f(P2) e ^ ^ ^ dSjdS2 


S S 


= livC^sH^ 


TT X .2 


^PO max ■ 1 


max 


If the antenna is assumed to be matched to the load. 


^o'^V^V max P*Pf 


RO 


2|Zi 


is the, power received on beam axis. Therefore, 


Pro 



ly 

2no 


2 

- = incident power density. 



AgQ is given in (A1.7). 


(Al.4.6) 


(Al.4.7) 


(Al.4.8) 


(Al.4.9) 


(Al.4.10) 
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If we expand 


® 2- m! ’ 

m=0 


(Al.4.11) 


then 


w_ = f ^ 

’po*^s^ m=0 




(Al.4.12) 


where 




Ipo(ks) 


(Al.4.13) 


Ipm(^s) = j J f(^i)f(P2) e 

S S 


. * -JV(W 


dSj^dS2 . 


(Al.4.14) 


Then 


-4 


2m 


9d = * 

m=0 


(Al.4.15) 


Icm(^s^ is called the correlation integral across the aperture. 
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A1.5 Signal Variance 


The signal variance is normalized to the measured dc power. 


Synchronous Receiver 


s 


2 

syn 






(Al.5.1) 



(Al.5.2) 


s^ dB = 10 log,« s^,,_ 
syn ^10 syn 


(Al.5.3) 


Using (Al.4.15) we may show that 


2 

^syn 



(Al.5.4) 


Asynchronous Receiver 


s 


2 

asyn 



(Al.5.5) 
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(Al.5.6) 


A1.6 


A1.7 


where 


9d 


asyn 


- 1 


4o^ 


- 29d + e X) 


Average Signal Level Degradation 



10 1 


P 

O9io p 


dc 

dc 


No fluctuations 


Synchronous Receiver 


AP, 


dcs 


4.343 o'" 


dB 


w 


Asynchronous Receiver 


AP, 


dca 


10 1 


dB 


O9io 


4ov 


|(1 + g<j)4(l - + e 


Average Effective Aperture 
The average effective aperture is 


- AgQgQ{k^)g^ , 


(Al.6.1) 


(Al.6.2) 


(Al.6.3) 


(Al.7.1) 


156 



(Al.7.2) 




I 


PO max 



(Al.7.3) 


AgQ is the maximum effective aperture for the plane wave case. 


A1.8 Average Directive Gain 

The average directive gain is 


Sd(l<s^ - Sq 9o(l<s^^d 


G 


0 max 





PO max * 


If 

Gq(^s^ = S^O^^s^ » 


then 

= S(^s)9d • 


(Al.8.1) 


(Al.8.2) 


(Al.8.3) 


(Al.8.4) 


(Al.8.5) 
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Gq niax maximum directive gain for plane wave incidence. Gq( 1<^) is 
the plane wave incidence power gain function, is the gain degradation 
factor. Finally, 


^d^’^s^dB “ ^0^’^s^dB ^°9l0 9d 


(Al.8.6) 


A1.9 Gain Degradation 

The gain degradation is defined as 

^^^^s^dB "" S^^s^dB ■ *^d^^s^dB * 
Therefore, 

^^(*^s^dB ~ ^°9l0 9d ' 

If (AI.4.15) is used. 





(Al.9.1) 


(Al.9.2) 


{Al.9.3) 
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A2.1 


From 


APPENDIX A2 

LOG AMPLITUDE AND PHASE CORRELATION 


Gaussian Spectrum 


The Gaussian Spectrum is, (3.3.3), 

2„2 

2 3 < ^n 

^n 


(3.5.2) - (3.5.4), in a transverse plane 


B |(p) = 2i\^k^L 


00 


2 3 po 

2. 2 On ^ 


2„2 

< Jin 


2TT^k"L 


0 


£ 

4 


/ 


0 


2 , 1 


_ sin(*^) 
1 ^ ^ 
kU 


2«2 
K Jlp 


6 Jq(kp)<CI< 
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The + sign is associated with B^(p) and the - sign with B (p). Separating 
the integrals, the first integral is. 





0 



from [ 30 ], ^ 6.631.4. 


The second integral is 


(A2.1.1) 


oc 

/ 


sin(^) 


2„2 

K Zn 


jQ(Kp)icdic . 


(A2.1.2) 


This is evaluated in Appendix Bl.l and is 



(A2.1.3) 


Combining (A2.1.1) and (A2.1.3), we have 



(A2.1.4) 

where 
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(A2.1.5) 


W = 


is called the Gaussian wave parameter. We also define the wave correla- 
tion function as 


B„(p) = B^(p) + B^(P). 


(A2.1.6) 


therefore. 


B„(o) 


= 


4 

2 , 7 , ^n 

TT e 

n n 


(A2.1.7) 


The Case of o = 0 


From (A2.1.5) and (A2.1.8), the variances are 






(A2.1.8a) 


(A2.1.8b) 


4 ' V” • 


(A2.1.8C) 


We note, finally, that if 'b' denotes the correlation coefficient, 
from (A2.1.6) and (A2.1.7), 




(A2.1.9) 
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A2.2 Von Karman Spectrum 


The Von Karman spectrum is 


= 0.033C^(4+ <) 


n 

6 



From (3.5.2) - (3.5.4) 



2 


Substituting t 



and defining 


2 

•^L 

2 

•^m 


_ 1 
2,2 
*^m^ 


W 


0 



» 


then 



(A2.2.1) 


(A2.2.2) 
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B ,(P) 

x»r 


o.osa/c^LQ^k^L 


00 

/ 

0 •- 


_ sin(WQt) 


(1 + t) « 


-U .7 t 




dt 


The integrals are evaluated separately. The first integral is 


“Znt r ^ r “Zgt a~l YTX“1 

/ e ^ (1+t) 6 Jg^'^LPAr^ dt = / e t (1+t) 


Jg(iCLP-y^) dt , (A2.2.3) 


_ 1 


where a = 1, y~ f • This is evaluated in Appendix (Bl.2.1) and is 


Z ” f-l)P /<|P\ 2p 1 

-^2 (-r) r(p+i)'f(p+i. P^. Zg) 

-n ( P* / 


p=0 


where y is a degenerate hypergeometric function. 


The second integral is 


f- 


sin(Wgt) -Zgt -11 

W^t ® Jq^'^LpV^^ • 
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Using sin(WQt) = Im 




we have 


1 I r "(Zn-jUn)* Y-a-1 J- 1 

rlmj / e “ t (1+t) Oo(<|_pVt) • 


where a = 0, y = -g- . Let 


^ " ^0 ■ • 


(A2.2.4) 


Again, using (Bl.2.1), 


^ J ^ r(p)'i'(p. P+'Y, z)l 

^p=o(P') ^ 


Applying Kummer's transformation [31] #13.1.29, 


1-b 

Y(a,b,Z) = Z f(l+a-b, 2-b, Z) 


and substituting -g- fory, the integral becomes 


1^ 2 


p^o ^ 


Im 


— -P 

r(p)z ® 'f(^, i^-p, z) 
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Therefore, 


B,(p) 

B*(p) 


5 

O.OSSiT^C^LQ^k^L 



/ 

< r(p+i)y(p+i, p-4’^0^ 


+ ^ Im 

”o 


11 


r(p)z ° ^-p, z)J 


where 


Z = Zn - jW, 


0 '^”0 


(A2.2.5) 



1 


S'-o 




(A2.2.6a) 


u Ho fl: 

~ Zq " k 


(A2.2.6b) 


Wq and may be defined as wave parameters for the Von Karman spectrum 
Next, using (A2.1.6), the wave correlation function is 


5 

B^(p) = 0.033(2 tt^)C^LqVl 



r(p+l)'f(p+l, p-t^, Zq). 
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For p = 0 


Only the p = 0 term exists in (A2.2.5). Then, applying (Bl.2.2) to 
the second degenerate hypergeometric function in (A2.2.5), we obtain 




0.033 ‘ 4'd. Zq) Im 






+ r(o) 


Now, since k^Lq » 1, Zq -► 0. Hence, from [31] #13.5.10, 


’’<^•'=>^0)* (0 < b < 1) 


so that 


'fd, 


1 

6 ’ 



6 

5 ’ 


(A2.2.7) 


Note that r(0) is pure real and, hence, drops out of the results. 
Performing the indicated substitutions and after some algebraic manipula 
tions, we obtain the final solution. 


5 

/ 

" .11 

= 0 . 033 / 1 ^ 

^ ^ ^ I'" 

(-1) ^y(^. 

1 __ 



i 


(A2.2.8) 
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(A2.2.9) 


al = 0 . 033 tt ^ * 

Once more, note that in general. 



A2.3 Kolmogorov Spectrum 


If the Kolmogorov spectrum (3.4.3), 




2 

JL- 

.U 2 

0.033C^ K 3 e 


were used instead, a solution for is not possible 

2 ^ 

is singular. However, a may be obtained since f (<) 

. . ^ X 

origin. 


o^ = 0.033C^(27T^)k^ J 


0 *- 


1 ■ 


-1 




u = 

U -3 . 


The first integral, from (Bl.3.3), is 


r« + 1) , 

— T— (4)' 


since the integrand 
is zero at the 


y+1 

d< 
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The second integral, from (Bl.1.7) is 


u -1 

.-Sin{ftan WJ 

2 WJl.w2)'> 


where 


W„ = 


<^L 

m 


m k 


is the wave parameter for the Kolmogorov spectrum. Thus, 


a 

/ 


0 •- 


1 - 


s1n(£ji) 


-• .>L- 

e K die 


r(-^ + 1) 


,P+2 

'm 


2 sin(^ 


1 - 


M 

P(l+w2)4 


•m 


and 


o = 


0.033ii2 I r(5)c2<'3k2L 


- , ji sin(ii tan'^W 

t! = 


m 


(A2.3.1) 

- 1 . 
(A2.3.2) 


168 



APPENDIX A3 


MOMENTS OF THE RECEIVED SIGNAL 


The assumptions implicit in this discussion are stated at the 
beginning of Section (3.8). The analysis presented in this appendix 
follows that of Shifrin [11]. 

Since it is assumed that the polarization vector of the incident 
wave does not vary with the position on the aperture, only the 
scalar components of the wave will be considered in the following. 


A3.1 First Moment of the Aperture Field 
As in (2.5.4), 

_ x(P)+J$(P) 

U(P) = Aji e 

■>. _ x(P)*j$(P) 

< U(P) > = Aq < e > . (A3. 1.1) 

% % 

Chernov [4] has shown that, for a Gaussian turbulence spectrum, x aid 
are jointly normal. We may evaluate (A3. 1.1) by the method of character- 
istic functions [32]., For a four-dimensional normal variable (xj^, X 2 » 

^ 3 » ^ 4 )* the characteristic function is 
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g(ti, 


^ 3 ’ ^ 4 ^ ^ 


< e 


j ( t^Xj^+t2X2+t2X2*^t^x^ ] 


4 1 4 

J I, 7 Z 

3 k=l 1=1 


4 

E Oj^aibkitj^ti 

k=l 


Let 


(A3. 1.2) 


Xj = X(P), X 2 = i|)(P). 


Set 


t 


1 " 


-j- 



Then 


< U{P) > = Aq g(-j, 1, 0, 0) 

Xo(F)*^[oJm*2Ja^(F)o^(F)bj^(F.F)-c,2(F)] 

Aq e 

where 

< X(P) > = XqCP) 

< 4)(P) > = 0 . (A3. 1.3) 


The medium is homogeneous. Thus, no point on the aperture plane 
has statistical privilege. Hence, 
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A 

V 

II 

(A3. 1.4a) 



(A3. 1.4b) 


Q 

II 

« 

(A3. 1.4c) 

Therefore, 



< U^(7) > = A„ . 


Let 

Uo = A„ . 

(A3. 1.5) 

Using 

(3.5.13b), we may then write 



U 

< U(F) > = Uq e ^ . 

(A3. 1.6) 

A3. 2 

Cross Moment of the Aperture Field 



< U(Fj) U(F 2 ) > = <e ^ ^ L i dj 

(A3. 2.1) 


Again, using characteristic functions with 
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= xCPj), X2 = X(P 2 ). X3 = J(P^). X4 = $(P2) , 

we find that [32] 

< U(Pi) U(P 2 )*> = g(-o, -j, 1. -1) 


Xo(Pi)+Xq(P2) - 

■ ^0 ® 


where 




2jax^cJ^^bx2’l^i + - 2a^^a^(,^b^|,^^^ 


The shorthand notation 1 and 2 for and P’ 2 , respectively, has 
been adopted for convenience here. Fortunately, under the assumptions 
made in Section (3.8), this reduces to 


L 





Therefore, 


u(v^) 


> = 


An G 


2Xo^X -[°x ' -’' 2 ) )] 


Using (3.5.13c), this reduces immediately to 
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< U(F^) U(F2)' 


> = An e 


2Xo*20x 


- 0 ^( 1 - 


"nC” 


I.P2)) 


Finally, from {A3. 1.5) 


< U(Pi) U(P2)* 


l“ol‘ 


-a 




w 


(A3. 2. 2) 
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APPENDIX A4 


ASYNCHRONOUS RECEIVER CHARACTERISTICS 


The fourth moment, dc power, and the variance of the signal volt- 
age measured by an asynchronous amplitude modulation (AM) receiver will 
be evaluated in this appendix. The receiver will be assumed to incorpor- 
ate a square law detector. It is, therefore, an envelope detector. 


A4.1 The Fourth Moment of the Measured Voltage 


We seek to evaluate the fourth moment of |V^(. Using (4.1.4) 


i4 ^ _ \u i4 


< ivr >= |l<vl IP'I’fl < / / e 


4 , /T ,2 , * 


//' 

s s 




•f's-C’!-'’?) , ff ,2 X3+VJ('>'3-*4) . » 

dSidSj J J ^ ^ 


s s 




e " ■ dS^dS^ > 


(A4.1.1) 
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where the notation is used to represent xCPj). f(P]^)» 

etc. If the statistical expectation operator and the integrals are 

interchanged, we obtain 


< 





/\ /\ y. y. 

iP-PfP ''o 


////< • 

s s s s 


Xi+X2+X3+X4+j('lti-4^+'l^-'l^ 4) 


V 2 V 4 ^ 


dS2dS2dS2dS^ 


(A4.1.2) 


Since Xj and i|ij are Gaussian random variables, we use the method of char- 
acteristic functions as in Appendix A3, 

8 

J t “ktk 

9^ > ^2 * ^3 ’ ^4 ’ ^5 ’ ^6 ’ ^7 ’ ^8 ^ - < 6 > 


8 

E 

k=l 




8 8 

E E 

k=l k=l 


^k'^l^kl Vl 


(A4.1.3) 


Let 


1^ ^ ^Xi+X2+X3+X4+j('l'i-'1^2+^3-'i'4) 


(A4.1.4) 


Assuming no statistical privilege, let 

< Xk > = Xo 


(A4.1.5a) 


2 

= a 

X 
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(A4.1.5b) 



<n> = 0 


(A4.1.5C) 



(A4.1.5d) 


We evaluate the function 


M 9 ( “ J > ” J » ~ J > ^ ^ ~ ^ ) • 


(A4.1.6) 


The first set of terms in the exponent is 


8 

^ ^ ^ ^ + <X2>+<X3> + <X4> 

k=l 

= 4 Xq . (A4.1.7) 


Let the second set of terms be written as 
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8 8 




k=l 1=1 


O a ^ ■*■ CTv V "*■ ^Y Y ■'' *^Y °Y ^Y Y 


* * ws^ 


+2 




^ ^\|j lb ~ *^\b ^ib ^4) ip 

^2 ^4 W . ^3 ^4 ^3 4 


+2j 


^ - \WS ^ \Xh\ 


'W\\\ * WhXz ' \X\h * \Xh% 


■'’X3“t';'’x34'i * '’X3‘’1'2‘’X3»2 ‘ W\*3 * %%'’!<3''’4 


°X4°*l'^X4*l * '’X4''*2'’x4't2 ' ‘’x4%'’x4>l^ * °X4%*’x4>l'4 


(A4.1.8) 
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The assumption has been made in the above that 


for any pair k,l. After simplifying and using (A4.1.5), this reduces to 
give 


* <4 


by -yj bv Y by V by V V 

^1 c. ^1 O ^4 C O t 


+ b 


X3X4J 


+ o* 


^ll) lb " bjj ^ + b,h (h ” lb 

^V2 ^r3 ^r4 ^2^3 ^2M 


^3*^ * 2J0v0.p. I >>•• ■. - ‘>. 


x"n >'i'i'3 ■ 


(A4.1.9) 


Combining the two sets of terms (A4.1.7) and (A4.1.9) and using the set 
of relationships (3.5.13), we get 


2 

Xi+X2+X3+X4+J(’J'i-'J'2+V'^4^ , 4Xq+4o^ ‘^w^^^'^n^ng^ 

^0 ^ 0 0 


"^w(l"^3n4^ °X^^iX3^‘\iX4^'\2X3''^2^^ 


^'"4'^‘''1'i'^4'‘'4^1'^3'‘^'1'24>3‘V4^ ^^'^‘"xV\i4^"‘'x2'1'4^ ^ (A4.1.10) 
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Clearly, the evaluation of (A4.1.2) using (A4.1.10) is impractical, if 

not impossible. Simplifying assumptions need to be made to obtain mean- 

4 

ingful results for <|V^ | > . 


We assume that the correlation of the refractive index is very 

high over the aperture of the antenna. Thus, if b^ tends to 1, then 

the log amplitude correlation coefficients b^ ^ become approximately 

X^Xi 

equal to each other, as do the phase correlation coefficients b, , . The 


cross correlation coefficients b , 

xi'i'a 


and b 


X2l>4 


also approach each other. 


Consequently, under these assumptions 


4Xo+4a^ Aal 
M = e X e X 


(A4.1.11) 


Using this in (A4.1.2) and using (3.8.1) gives 



\%\ 1^1 |P*Pfl 


•••'■ ////v 

s s s s 




f 3 f 4 e ^ ^ dS^dS 2 dS 3 dS 4 


which may be reduced, by using the definitions in Sections (4.2) and 
(4.3) to give 


< iVp 1^ > = iVjol^ |gvC<s) I"* |P*Pf I 


4o? 


17 ^ In.n e ^ 


(A4.1.12) 


179 



Finally, it is easily shown, using (A4.1.12) that 



7 — ? 4 

' ’’rO 3of's> IP-M ® 


(A4.1.13) 


A4.2 Signal Level 


The dc power measured by an asynchronous receiver assuming a match- 
ed antenna load is 


^dca 


2|Z,l^ 



(A4.2.1) 


< |VJ 


2|Zi 


T - '^1 



Rl > 


2 


% 

The square law detector precludes the measurement of |V^ | directly. 
Also, the total power is (4.3.9), 


“ *^R0 Sq^'^s^ |P‘Pf 1^ 


' 9d » 


(A4.2.2) 
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where 


^0 " '’rO 


(A4.2.3) 


The arguments have been dropped for convenience. Then 


Pdca ^ 


For small fluctuations, 



From Taylor's theorem, for x near 1, 





= -^(1+x) - •^(l-2x+x^) X * 1 


(A4.2.4) 


Therefore, 
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< 





1 - 


2 > 


< ijlM 



(A4.2.5) 


Squaring both sides and using (A4.1.13), we obtain 


^dca ~ ^RO^O^'^s^ 


/ 

|(l+9d) 


4a'' 

|(l-2Ve X) 


^ /s. O 

|P*Pf| 


(A4.2.6) 


Finally, taking the square root and simplifying this gives 


asyn 


l''dol |9vC‘s>l 


y 2 ' 

■ |(i*9d) - i(i-2g/e'''’x) 






lp.pfl* 

(A4.2.7) 


A4.3 Variance 


From (4.1.8), the normalized variance measured by an asynchronous 
receiver may be written as 


s 


2 

asyn 


■ ^dca _ ^R 

P ~ P 

dca dca 


(A4.3.1) 
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Substituting from (4.3.8) and (A4.2.6), we obtain 


s 


2 

asyn 



fd , 

^(l-29d+e X) 


- 1 . 


(A4.3.2) 
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APPENDIX A5 

EVALUATION OF THE CIRCULAR APERTURE 


Uniform Illumination 

The uniform aperture illumination function: f (R) '= 1. 


A5.1 The Voltage Gain Function gy(0) 


The integral ly, from (4.2.3) » becomes 



1 2n 


// 

0 0 


j0RCOS(C-(j)') 5 

e Rd^dR = 


I 


Jq( 0R) RdR 


From [30] #6.561.5, 



\H-1 

X 


0^(ax)dx = 


Vl(a) 

a 
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Therefore, 


2Ji(0) 

Iy(®) = %hys ~ ’ 


(A5.1.1) 


where 



is the physical area of the aperture. 


I 


V max 


= A 


phys 


(A5.1.2) 


Hence, 


2Ji(0) 

9vO) = -T- 


(A5.1.3) 


The normalized power gain pattern function is 


90 <®> 


2Ji(e) 

~0 


(A5.1.4) 
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A5.2 The Integral Ij,(C,e) 


From (4.3.10) - (4.3.13), 


I(C, 0 ,(j>) = 


// 

S S 


— 2 ^ 

_ C 


jeR^cos(Ci-<t>) 


-j0RpCos(5p-(j>) 

e RjR2dqdR2d^dR2 


Using the power series expansion for the exponential 


KC.e.tj,) = a^ 


p=o\^/ ss 


R^+R 


2 

2 



j0R,cos(^,-<|)) -j 0 RpCOs(^p- 4 )) 

e ^ ^ e R^R 2 d 5 ^dR^dC 2 dR 2 . 


Now 


“^”<51-52) (0 


cos 


(p-2q)(q-?2) 


q=0 
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Considering the ^ integrals only, 


H TT 


■//[ 

“IT “TT 


cos(p-2q)q cos(p-2q)^2 * sin(p-2q)q sin(p-2q)^2 


j0R,cos( q-(j)) . -j0R2COS(^2"4>) 
e e dqd^2 « 


but 


/ 


sin(p-2q)5 


+j0Rcos(5-(|)) 

e 


dc 


-n 


0 


so that 


/ 


Q = I cos(p-2q)Cj e 


j0Rj^cos(5j-(j,) 


dCi 


-TT 


/ 


-TT 


cos(p-2q)C5 


-j0R2COs(52"(l)) 
e ^^2 * 


From [30] #3.915.2 


/ 


jbcosx 

e cosnxdx = Tt J^(b) 


> 
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therefore. 


(-l)” %n^Jp.2,(0Rl)Op.2q(-0R2) 


but 

Jp(-x) = (-1)%(X) , 


thus 


and is independent of 4>. Then 


KC.0) 





0 


Hence, 


/, 
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(A5.2.1) 



Define 


f "c^ P+1 

Qp = j e ^ R Jp_2q(0R)dR • 
0 


(A5.2.2) 


From Appendix (B2.1.2), using 


V = p-2q 


(A5.2.3) 


]i = 2t+p+l 


(A5.2.4) 


Let 


P+2 

Qp = V" 


and 


6 


C 


C0 


(A5.2.5) 


(A5.2.6) 
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where 


Qc(C,0) = 


(- 1 )^ 1 


Z! t! g^vtfl 
t=0 X 


ej (0)s , , 

2U V u-l.v-1 


^ r(^)J 


(A5.2.6) 


Therefore, 



(A5.2.7) 


A5.3 On-Axis Incidence 

For the special case of 0 = 0, 

9v(0) = 1 

9o(0) = 1 


(A5.3.1) 

(A5.3.2) 


Iq(C, 0) is non-zero only if Qq(C, 0) is non-zero. From (A5.2.1) and 
(A5.2.6) with 


p-2q = 0 


Qc(C.O) = Y(q+1. . 
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APPENDIX A6 


EVALUATION OF THE CIRCULAR APERTURE 


Gaussian Illumination 


The Gaussian aperture illumination function: 


f(^) = 



A6.1 The Voltage Gain Function 9 ^( 0 ) 


The integral ly, from (4.2.3) becomes 


I 


V 



■_2 j©Rcos(C-<t>) 
e ^ e Rd?dr 


1 



\2 

e JQ( 0 R)RdR . 


Using the results of Appendix B2.1, substituting x for C, n=v=0 and 


T0 
2 . 
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(A6.1.1b) 


2 

"" ^^phys f" 







Using (Bl.3.2), when 0 = 0, 


max 








/ 

r 

1 

^PO max 


1 

' Q} 

1 





(A6.1.2a) 


(A6.1.3a) 


(A6.1.4a) 
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The Case of Small t 


If the aperture taper constant t is small, so that the illumination 
beyond the edge of the aperture is small, the integrals may be taken 
to infinity with small error. Then, using [30] #6.631.4 


’vl®) = \hys-'^ ® ^ 


(A6.1.1b) 


max “ %hys'^ 


(A6.1.2b) 


2^2 

_ T a 

0^ 4 

gv(e) = e ^ = e 


(A6.1.3b) 


A6.2 The Integral Iv(C,0) 


From (4.3.12), let 


I(C,T,0,(J>) “ a 




r 2 ,^2 

e ^ 


S S 


j0RiCos(5i-(t>) -j0RoCOs(^p-())) 

e ^ ^ e ^ RjR 2 dR^dR 2 dCid ?2 • 


Defining 



1 

T 

y 
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and performing the ^ integrals gives 


((C.T.e) . 



(A6.2.1) 



/ 


'y2 p+1 


(A6.2.2) 


Again, using (82. 1.2), and the results of A5.2, 


V = p-2q 


(A6.2.3) 


U = 2t+p+l 


(A6.2.4) 


Let 



vP+2 

= 2 “ %(y^Q) 


(A6.2.5) 



(A6.2.6) 


where Q^(y,0) 
\for 


is the identical to (A5.2.6) 


with Y substituted for C and 
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Hence, 





(A6.2.7) 


A6.3 On-Axis Incidence 

/ 

For the special case 0 = 0, 


9v(0) = 1 


(A6.3.1) 


9o{0) = 1 


(A6.3.2) 


I (C,0) is non-zero only if Qp(y»9) f 0. From (A6.2.2) we have 
Y ^ 


p - 2q = 0 


Qc(Y,0) = Y(q+1, . 
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APPENDIX B1 


EVALUATION OF INTEGRALS 


Bl.l 


CO 

■/ 


sin(AK^) 


A<^ 


J^(p<) d< . 


(Bl.1.1) 


Using 


sin(Aic'^) = Jm 


' JAK^' 

e 

V t 


Let 


C = 4(B-jA) 


(Bl.l. 2) 


Then, 


I 



JJpk) d<' 
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The integral can be evaluated directly using [30] #6.631.5 to give 


I 


"sMAEii 


A< 


2 e < J^(P<) dK 


' 2Ar(v+iy "gu * iV • '<+1. -r j 
I c ^ J 


where $ (a,b,Z) is a degenerate hypergeometric function. 


The Case of u=u 


From [31] #13.6.10, 


$(a, a+1, -x) = ax”S'(a,x) . 


Using this in (Bl.1.3), we get 


I p = V 


Im 


/ 


k-jk+a 


^ (-irr 

where yict,Z) = ^ ' klC^y 


k=0 


is an incomplete gamma function. Therefore, 


(Bl.1.3) 


(Bl.1.4) 
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The Case of u=v=0 



Also, in (Bl.1.4) [30] # 8.356.3 and # 8.359.1, 


y(a,Z) = r(a) - r(a,Z) 


r(0,Z) = -Ei(-Z) . 


Therefore, 



f 

1 

o\ 




ly=V=0 


> 


Where r(a,x) is an incomplete gamma function 
Ei(x) is the exponential integral. 


The Case of v=0. P=0 
For this case (61.1,3) reduces to 


f ^ 



(61. 1.6) 


(61.1.7) 
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Thus, after some algebraic manipulations, 


/ 


sin(A<^) 



-Btc^ M+1 
e ic die 



(Bl.1.8) 


B1.2 

Consider the integral 


f -It a-1 Y-a-1 ^ 

I = / e t (1+t) dt . 


Using the series expansion for Jq(x), we have 



-Zt p+a-1 p+Y-(p-Kx)-l 
e t (1+t) 


dt . 


From C 30 ], #9.211.4, this gives 


I 


y . -^^^2 r(p+a)'l'(p+a, P+Y» Z) . 

p=o 


(Bl.2.1) 
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where [30] #9.210.2, 

f [a-b+l) +li^Z^ ^(a-b+1, 2-b, Z) (Bl.2.2) 

and $(a,b,Z) and ’i'(a.b.Z) are degenerate hypergeometric functions. 

B1.3 

oo 

f -px^ 21+1 
I = j e X dx . 

u 

2 

Substituting y = x . we have 

OO 

I f e-^y y' 6y. 

From [ 30 ] #3.381.3, this gives 



-px^ 21+1 
e x dx 


u 


2p‘ ^ 


Similarly, using [30] #3.381.1 and #3.381.4, 


(Bl.3.1) 
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i,:j-T(l+l, pu^) 


(Bl.3.2) 


/ 


-px^ 21+1 
e X dx = 


2p 


I 


-px^ 21+1 
e X dx 



r(i+i) 


(Bl.3.3) 


y(a,x), r(a,x) are the incomplete gamma functions of the first and second 
kind, respectively. r(x) is the gamma function. 
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APPENDIX B2 


B2.1.1 


To evaluate 



From [30] #6.561.13 
1 

J x^J^(0x)dx 

0 


EVALUATION OF INTEGRALS 


2n+v+l 

J^(6R) dR 


(B2.1.1) 



2{m+n)+v+l 

R J^(0R) dR . 


1 

0P+1 




(e)-0J i(0)s (0) 

\rt y»v 



Re(\r+y) > -1 • 
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S (0) is a Lommel function. 
r(x) is the Gamma function. 


Note that the above result is the corrected form of the expression given 
in [30] ^*6. 561. 13. 

Here y - 2m-<-2>nv+l. Multiply and divide by Then, 


Qn = C 


2n+v+2 


f* (- 1 )" ( 1 Y*' 

L n,. [eg; 

m=o 




r(^)l 

J,^,(e)s (0) + 2*' 

v-1' y,v' v-u+1 


R2.( vfy) > -1 . 


(B2.1.2) 


B2.2 For 0 = 0 


From B2.1.1, Qp can now exist only if v = 0. Therefore, 


1 


Qp(0) “ y* 


■p2 2n+l 
e '' R dR. 


0 


From Bl.3.2, this gives 


p2n+2 1 

Qp(0) = Y(n+1, . 

C 


(B2.2.1) 
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